!!'*| | l , !"p. ! ! | 1 1 1 -j! !'=» pji'l'ff ; t!ir^[r i?, f ITI'fl'f] '^1 r "* fl 




NASA TECHNICAL 

MEMORANDU 


REAL-TIME SIMULATI 
THE TF30-P-3 TURBG 
USING A HYBRID C 




by John R. Szuch and William 

Lewis Research Center 
Cleveland, Ohio 44135 


NATIONAL AERONAUTICS AND SPACE AD 




1 Report No. 2. Government Accession No, 

NASA TM X-3106 

3. Recipient's Catalog No. 

4 . Title and Subtitle 

REAL-TIME SIMULATION OF THE TF30-P-3 TURBOFAN 
ENGINE USING A HYBRID COMPUTER 

5. Report Date 

OCTOBER 1974 

6. Performing Organization Code 

7. Author(s) 

John R. Szuch and William M. Bruton 

8. Performing Organization Report No, 

E-7904 

10. Work Unit No. 
501-24 

9. Performing Organization Name and Address 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 44135 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Memorandum 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 


A real-time, hybrid-computer simulation of the TF30-P-3 turbofan engine has been developed. 
The simulation is primarily analog in nature but uses the digital portion of the hybrid computer 
to perform bivariate function generation associated with the performance of the engine’s rotating 
components. FORTRAN listings and analog patching diagrams are provided in the report. The 
hybrid simulation was controlled by a digital computer programmed to simulate the engine's 
standard hydromechanical control. Both steady-state and dynamic data obtained from the 
digitally controlled engine simulation are presented. Hybrid simulation data are compared with 
data obtained from a digital simulation provided by the engine manufacturer. The comparisons 
indicate that the real-time hybrid simulation adequately matches the baseline digital simulation. 


17. Key Words (Suggested by Author(s)) 

Simulation Real-time 

Hybrid computer Transient 

Turbofan Dynamics 

18. Distribution Statement 

Unclassified - unlimited 
Category 28 


19. Security Classif. (of this report) 
Unclassified 

20. Security Classif. (of this page) 
Unclassified 

21. No. of Pages 

103 

22. Price" 

$4. 50 


* For sale by the National Technical Information Service, Springfield, Virginia 22151 




CONTENTS 


Page 

SUMMARY 1 

INTRODUCTION 1 

ENGINE DESCRIPTION 3 

ENGINE MODEL 4 

Inlet 4 

Fan 5 

Low-Pressure Compressor 6 

High-Pressure Compressor 7 

High- and Low-Pressure Turbines 9 

Combustor, Augmentor, and Duct 10 

Exhaust Nozzle 12 

Engine Dynamics 13 

HYBRID SIMULATION i8 

Digital Program I 8 

MAIN digital program I 8 

MAP2 function routine 19 

Analog Program 19 

Digital-to-analog multiplier interface 20 

Analog function generation 20 

RESULTS AND DISCUSSION 20 

Steady State 20 

Transients 22 

Frequency Responses 24 

CONCLUDING REMARKS 25 

APPENDIXES 

A - SYMBOLS 27 

B - MAIN DIGITAL PROGRAM - DESCRIPTION AND FORTRAN LISTING .... 35 

C - DIGITAL INPUT DATA 41 

D - MAP2 FUNCTION GENERATION ROUTINE - DESCRIPTION AND FORTRAN 

AND ASSEMBLY LANGUAGE LISTINGS 43h 

E - ANALOG PROGRAM - PATCHING DIAGRAMS, POTENTIOMETER SETTINGS, 

AND SCALE FACTORS 51 

REFERENCES 52 

iii 



REAL-TIME SIMULATION OF THE TF30-P-3 TURBOFAN ENGINE 
USING A HYBRID COMPUTER 
by John R. Szuch and William M. Bruton 
Lewis Research Center 

SUMMARY 

This report describes the development of a real-time hybrid-computer simulation 
of the Pratt & Whitney TF30-P-3 augmented turbofan engine. The simulation is intended 
to support research programs involving this engine. 

The simulation is primarily analog in nature but does use the digital portion of the 
hybrid computer to perform bivariate function generation associated with modeling the 
performance of the engine's rotating components. The digital portion of the engine sim- 
ulation requires about 10 000 words of core storage. The digital sampling interval is 
about 4. 6 milliseconds and does allow real-time simulation of engine dynamics. Ap- 
proximately 185 amplifiers and 85 multipliers are required to perform the analog calcu- 
lations. FORTRAN listings and analog patching diagrams are provided in appendixes to 

this report. 

The hybrid simulation was controlled by a general-purpose digital computer pro 
grammed to simulate the standard TF30-P-3 hydromechanical fuel and variable geom- 
etry controls. Both steady- state and dynamic data obtained from the digitally controlled 
engine simulation are presented. Hybrid simulation data are compared with data ob- 
tained from a digital simulation provided by the engine manufacturer. These compari- 
sons indicate that the real-time hybrid simulation does adequately match the baseline 
digital simulation. The real-time simulation should, therefore, provide a valuable 
analytical tool for research studies involving the TF30-P-3 turbofan engine. 


INTRODUCTION 

Aircraft propulsion systems have become steadily more complex. With engines and 
inlets operating at higher levels of performance, propulsion system controls have also 
necessarily increased in complexity. In order to obtain maximum performance for all 



flight conditions, digital control systems (ref. 1) are being considered for inlets and 
engines. The use of a digital computer, with its logic and memory capabilities, would 
provide maximum flexibility and versatility to propulsion control systems. 

The development of controls for aircraft propulsion systems is facilitated by an 
ability to predict engine performance accurately over a wide range of operation. Com- 
puter programs for analyzing steady-state and dynamic performance of turbojet and 
turbofan engines have previously been developed (refs. 2 to 8). Some of these engine 
simulations have been expanded to include the simulation of engine controls. 

In the development of digital controls for turbojet and turbofan engines, both com- 
puter simulations and experimental testing will be required to evaluate new control laws. 
The experimental tests can be more efficiently planned and executed if the digital control 
software can be tested and ’'debugged'’ by using the engine simulation. This imposes an- 
other requirement on the engine simulation. That is, the simulation must operate in 
"real time. " The real-time requirement exists because timing and sequencing are ex- 
tremely important in the development of digital control software. Also, the effects of 
digital sampling intervals on engine performance (ref. 9) can be predicted by using a 
real-time simulation. 

Turbojet and turbofan engines can be simulated by using either analog (ref. 6), 
digital (refs. 2 to 5), or hybrid (refs. 7 and 8) computers. The real-time simulation 
requirement, however, limits the possible choices to the analog and hybrid computers. 

A purely analog simulation of a turbofan engine requires a large amount of computing 
equipment. Much of this equipment is needed to perform the bivariate function genera- 
tion associated with modeling the engine’s fan, compressor (or compressors), and tur- 
bine (or turbines). Unfortunately, analog function generation is often plagued by a lack 
of precision and repeatability. The precision and memory capabilities of the hybrid’s 
own digital computer can be used to perform the required function generation, thus re- 
ducing the required amount of analog computing equipment. 

This report describes the development of a real-time hybrid-computer simulation 
of the Pratt & Whitney TF30-P-3 augmented turbofan engine. The simulation is intended 
to support research programs involving this engine. The hybrid simulation was pat- 
terned after the engine manufacturer's digital simulation of the TF30-P-3 engine and 
was implemented on the Lewis Research Center’s Electronic Associates, Inc. , (EAI) 
Model 690 Hybrid Computing System and two EAI Model 231-R Analog Computers. 

The TF30-P-3 hybrid-computer simulation was controlled by a general-purpose 
digital control computer (ref. 10) programmed to simulate the standard TF30-P-3 
hydromechanical fuel and variable geometry controls. This digital control system was 
developed for on-line control of airbreathing engines during experimental testing. 

Both steady- state and dynamic data obtained from the digitally controlled simulation 
are presented for a full range of power settings. Data were obtained at both sea level 
and simulated altitude conditions. Hybrid simulation data are compared with data ob- 
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tained from the engine manufacturer's digital simulation. FORTRAN listing and analog 
patching diagrams are provided. This information should prove useful in the develop- 
ment of similar engine simulations. 


ENGINE DESCRIPTION 

The Pratt & Whitney TF30-P-3 engine (fig. 1) is an axial, mixed-flow, augmented, 
twin-spool, low-bypass-ratio turbofan. A single inlet is used for both the fan airflow 
and the engine core airflow. Airflow leaving the fan is separated into two flow streams: 
one stream passing through the engine core, and the other stream passing through the 
annular fan duct. A nine-stage low-pressure compressor includes the core portion of 
the three-stage fan and is connected by a through-shaft to the three-stage low-pressure 
turbine. A seven-stage high-pressure compressor is connected by a hollow shaft to the 
single-stage high-pressure turbine. A seventh-stage compressor bleed system dis- 
charges air into the fan duct at flight Mach numbers above 1. 75. A 12th- stage bleed 
automatically discharges air into the fan duct for compressor surge control. Cooling of 
high-pressure and low-pressure turbine blades is accomplished by using 16th-stage 
compressor bleed air. 

The TF30-P-3 combustor is located between the discharge of the high-pressure 
compressor and the inlet of the high-pressure turbine. The combustor section consists 
of an annular diffuser and eight combustion chambers. Each combustion chamber con- 
tains four dual-orifice fuel spray nozzles. The main fuel control meters fuel to the 
combustor as a function of the power lever angle PLA, the high-pressure-compressor 
speed N h , the fan inlet total pressure P 2 and temperature T 2 , and the high-pressure- 
compressor discharge static pressure P g 3 . (All symbols are defined in appendix A. ) 

In the combustor, airflow from the compressors reacts with the injected fuel, producing 
high-energy gas to run the turbines. 

The engine core and fan duct flow streams combine in an augmentor and are dis- 
charged through a variable-area convergent nozzle. The augmentor consists of a dif- 
fuser section and a combustion chamber. The fuel is introduced into the combustion 
chamber through spray rings which are arranged into five separate zones. The augmen- 
tor fuel control meters fuel to the five zones as a function of PLA, N H , T 2 , P g 3 , the 
nozzle area A N , and the turbine discharge total pressure P g . The hydraulically ac- 
tuated nozzle is adjusted to maintain a scheduled ratio of P g 3 /Pg during augmentor 
operation. 
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ENGINE MODEL 


The first step in developing an engine simulation is the formulation of an analytical 
model. This model, in equation form, represents the functional relations that exist be- 
tween engine variables such as pressures, temperatures, and flows. In the case of the 
TF30-P-3 hybrid-computer simulation, the engine model should be capable of accurately 
predicting both steady-state and dynamic performance of the engine. Wide-range per- 
formance maps for the engine's rotating components should be included in the simulation 
for steady- state accuracy. Factors such as fluid momentum, mass and energy storage, 
and rotor inertias should be included to provide the desired transient capability. 

Wherever possible, the TF30-P-3 model was patterned after the engine manufac- 
turer s baseline digital simulation of the TF30-P-3 engine. The following sections de- 
scribe the various elements which comprise the TF30-P-3 hybrid-computer simulation. 
Figure 2 illustrates the flow of calculations in the hybrid simulation and should help the 
reader in understanding the relations of the individual elements in the engine model. 


Inlet 

The real-time, hybrid- computer simulation of the TF30-P-3 engine does not include 
a dynamic simulation of an inlet. However, a steady-state representation of a typical 
inlet recovery is included in the simulation. The following equations are used in the hy- 
brid simulation to calculate the total pressure P 2 and total temperature T 2 at the fan 
inlet from the values of Pq, Tq, and Mq for a specified flight condition: 


(T/T)j = 1 + 0. 2 Mq 

(1) 

(P/P)j = (T/T)j ■ 5 

(2) 

a. 0 if Mq as 1. 0 ■> 


ll.O - 0.075(M 0 - l) 1 ’ 35 if M q >1. (L 

r (3) 

t 2 = (T/TJjTq 

(4) 

P 2 = (P/P)^ 

(5) 
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Fan 


There are basically two methods that can be used to model multistage fans and com- 
pressors. One method is to represent multistage fans and compressors with individual 
stage models (i. e. , to compute pressure and temperature rises across each stage). 

This technique is referred to as stage stacking (ref. 6), but it requires a large comput- 
ing facility when used in a total engine simulation. The method used in the real-time, 
hybrid-computer simulation involves the use of overall component performance maps. 
Interstage gas dynamics were not considered. 

The TF30-P-3 three-stage fan was modeled with separate performance maps for the 
fan hub (core) and tip (duct) sections. Separate maps were required because of the ra- 
dial pressure gradient which exists at the fan discharge at high rotor speeds. 

Overall fan performance for the hub and tip sections is described by a pair of bi 
variate functions. In the baseline digital simulation, the fan tip pressure ratio and fan 
corrected speed were selected as independent variables and used to compute the total 
fan corrected airflow and the fan tip adiabatic efficiency. The same pair of independent 
variables were used to compute the fan hub pressure ratio and fan hub adiabatic effi- 
ciency. 

To obtain the fan discharge temperature from efficiency requires additional calcu- 
lations (including exponentiation). Because of the real-time requirement and equipment 
limitations, it was necessary to minimize the function generation and algebraic compu- 
tations associated with modeling the fan (and compressors). Data obtained from the 
baseline simulation indicated that the fan hub and tip efficiency maps could be eliminated 
and the corresponding temperature ratios calculated as linear functions of the pressure 
ratios. In the TF30-P-3 hybrid-computer simulation, the fan model is described by the 
following equations: 


(w c ) 
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'P» 1- 3 


^P, 1-3 


w BL12 T 3 


(T/T) F jj^ - 0. 7588 + 0. 2412(P/P)pjjj 


(ID 

( 12 ) 


T 2. 1 T 2^ T//T ^FID (13) 

Figure 3 contains plots of the bivariate functions representing the performance of the fan 
(eqs. (6) and (8)). The pressure ratios, corrected speed, and corrected airflow have 
been normalized by their design values. Equation (11) indicates that the 7th- and 12th- 
stage bleeds are assumed to enter the fan duct immediately behind the fan. The primed 
station designation (e. g. , 1. 3') refers to the entrance to the mixing volume. 


Low-Pressure Compressor 

Overall performance data are used to describe the TF30-P-3 low-pressure com- 
pressor (4th to 9th stages). The compressor corrected airflow and corrected speed 
were selected as the independent variables and were used to compute the compressor 
pressure ratio. As in the case of the fan, compressor temperature ratio was calculated 
from a linear function of pressure ratio, thus eliminating the need for generating the bi- 
variate efficiency function. In order to better match the temperature ratio - pressure 
ratio relation over the range of operation, two linear functions were used for different 
ranges of pressure ratio. 

The modeling of the low-pressure compressor was complicated by the effect of the 
7th -stage bleed on the overall compressor performance. Based on information obtained 
from the digital simulation, the effect of the bleed on performance was represented by a 
shift in the compressor corrected speed used in the pressure ratio determination. The 
effect of bleed on efficiency (and hence, temperature ratio) was neglected in both the 
baseline digital and real-time hybrid simulations. 

In order to compute the low-pressure-compressor corrected airflow with the 
7th -stage bleed open, the following procedure was used: First, the corrected speed 
shift for a fully open bleed was determined from a function of the actual corrected speed. 
Linear interpolation was then used to determine the speed shift for a partially open 
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bleed. Finally, the pressure ratio was calculated from a bivariate function of the cor- 
rected airflow and shifted corrected speed. 

The 7th-stage bleed flow was assumed to be proportional to the low-pressure- 
compressor inlet airflow and was assumed to leave the core stream at the compressor 
exit (station 2.2). 

The TF30-P-3 low-pressure-compressor model is described by the following equa- 
tions in the hybrid-computer simulation: 


(P/P>L C = *3 


s. i 


I' ( N c) 


2 . 1 

P 2. 2' = P 2. l( P//p ) LC 


LC, M 


N I 

- - + XBL7 fn 

X 

'LC,M y0 2 1 

_V 0 2. 1 


(T/T) 


LC 


( N c) 


fO. 84 + 0. 1989(P/P) LC if (P/P) LC > 1. 16 "| 


,0. 515 + 0. 4689(P/P) lc if (P/P) LC ^ 1- 16 J 


(14) 


(15) 


(16) 


(17) 


2 . 2 


= T 


2 . 1 


(T/T) 


LC 


(18) 


= 0. 06 XBL7 Wg i 


(19) 


Figure 4(a) contains a plot of the bivariate function representing the low-pressure- 
compressor performance (eq. (14)) with the 7th-stage bleed closed. Figure 4(b) contains 
a plot of the corrected speed shift for the 7th-stage bleed fully open. The low-pressure- 
compressor map variables have been normalized by their design- point values. 


High-Pressure Compressor 

The mathematical treatment of the high-pressure compressor (10th to 16th stages) 
was quite similar to that of the low-pressure compressor. The basic difference was in 
the treatment of the 12th- stage bleed and its effect on the overall compressor perform- 
ance. Based on information obtained from the digital simulation, the effect of the bleed 
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was represented by a shift in the compressor corrected flow used in the pressure ratio 
determination. The 12th- stage bleed flow was assumed to leave the core stream at the 
compressor exit (station 3). The cooling bleeds for the high- and low-pressure turbines 
are also extracted at the high-pressure-compressor exit. These bleed flows were as- 
sumed to be proportional to the high-pressure-compressor inlet airflow. 

The high-pressure-compressor temperature ratio was calculated from a linear func- 
tion of pressure ratio. Data from the digital simulation indicated that the slope and in- 
tercept of the function were sensitive to the bleed. Therefore, the slope and intercept 
were assumed to be linear functions of the percent of fully open bleed. The TF30-P-3 

high-pressure-compressor model is described by the following equations in the hybrid 
simulation: 


( p /P)jjc - f 4 


(*<=)■ 


N 


hc,m’ 


H 


2.2 


P 3 ,= P 2. 2 (P/P)hc 


("c) 


HC, M 


W 2.2V^2 

6 2.2 


- XBL12 f 


10 


N. 


H 


V 0 2. 2 


(T/T) hc - (1- 0322 - 0. 185 XBL12) + (0. 12588 + 0. 08032 XBL12)(P/P) 


HC 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


T 3' = T 2. 2 ( T / T ) H C (24) 

BL12 = °- 055 XBH2 w 2 2 (25) 

W BLHT = °- 04744 w 2 2 (26) 


W BLLT = °- 01985 w 2 2 (27) 

Figure 5(a) contains a plot of the normalized bivariate function representing the high- 
pressure-compressor performance (eq. (20)). This plot represents the closed-bleed 
performance of the high-pressure compressor (XBL12 - 0). Figure 5(b) contains a plot 
of the normalized corrected flow shift for the 12th-stage bleed fully open. 
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High- and Low-Pressure Turbines 


In the baseline digital simulation, overall turbine performance is represented by a 
set of bivariate functions. For each turbine, pressure ratio and a speed parameter 
N/\/t were used as independent variables to determine the turbine flow parameter 
w-^T/P and enthalpy drop parameter Ah/T. Based on previous simulation experience, 
the digital program map outputs were normalized by the speed parameter for use in the 
hybrid simulation models. The TF30-P-3 turbine models are described by the following 
equations in the hybrid simulation: 
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Figures 6 and 7 contain plots of the normalized bivariate functions representing the high- 
pressure- and low-pressure-turbine performance (eqs. (28), (30), (34), and (35)). The 
constants K fiWHT and Kg WLT appearing in equations (32) and (37) denote the frac- 
tions of the cooling bleed flows that are assumed to do turbine work. 


Combustor, Augmentor, and Duct 

The performance of the combustor can be represented by a total pressure loss due 
to friction and a temperature rise resulting from the heat release from the burned fuel. 
The ratio of the total pressure loss to the combustor inlet pressure P„ was assumed to 

u 

be proportional to the square of the compressible flow parameter w^T/P at the com- 
bustor inlet. The following equation is used in the hybrid simulation to calculate the re- 
sulting combustor exit pressure: 


Pa- = P, 


K B w 3 T 3 


= P 3- 


K B W 3 V 3 

R A W 3 


(38) 


The use of the stored mass in equation (38) results in a savings in computing equip- 
ment. 

The digital program makes use of a number of bivariate functions to describe the 
combustor efficiency and temperature rise. Data from that simulation indicated that a 
simple energy balance assuming constant combustion efficiency could adequately de- 
scribe the steady-state temperature rise in the combustor. The following equation is 
used in the hybrid simulation to describe that process: 


(wT) 4 , 


_ w 3 t 3 + ^ 
C P,4 


HVF 


c 


P, 4 



( 39 ) 


The specific heats c p , 3 and c 4 , together with the efficiency r? B and fuel heating 
value HVF, were assumed to be constant. 
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The TF30-P-3 augmentor was represented by two sections having separate flow 
streams. In the baseline digital program, a static pressure balance between the 
streams was assumed at the entrance to the augmentor. However, only total pressures 
are used in the hybrid simulation. Data from the digital simulation indicated that the 
total pressure ratio Pg/Pj g was approximately constant along selected engine oper- 
ating lines. Therefore, in the hybrid simulation, the total pressure Pg and Pg are 
calculated from 


P 5 " P 6 “ K PR6 P 1. 6 


(40) 


Each section of the augmentor was treated in the same way as the combustor. The fuel 
flow to the core augmentor section was supplied from two zones, and the duct augmentor 
section was fed from the remaining three zones. The following equations described the 
TF30-P-3 augmentor model: 


p 7 - p 6 - k c ab^ w 4. 1 + W BLLT^ 2 ^ 

0 

(41) 

(wT)g = (wT)g 

(42) 

(wT)n - ^ (wT)g + HVF w F7 

(43) 

C P,7 C p,7 


(wT) 7 

t 7 = - : : — 

(44) 

w 4. 1 + W BLLT + W F7 


„ „ K DAB w 1.6 V 1.6 

P 1.7’- P 1.6 RW 

k A w 1. 6 

(45) 


(wT) 


1.7 




"P, 1-7 


1. 6 x 1. 6 


+ hvf w 
C p,1.7 


FI. 7 


(46) 


As in the combustor, the specific heats and combustion efficiencies in the augmentor 
were assumed to be constant. 

The total pressure drop in the bypass duct was modeled in the same manner as for 
the combustor sections. However, no temperature rise (or drop) was assumed in the 
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duct. The following equations describe the hybrid simulation duct model: 


P l. 6’ " P l. 3 


K P W 1. 3 V 1. 3 
R A W 1.3 


( 47 ) 


T l. 6' “ T l. 3 


(48) 


Exhaust Nozzle 

In order to match the analytical model used in the digital simulation, the TF30-P-3 
exhaust nozzle was modeled as two separate nozzles fed by the separate core and duct 
augmentor sections. The sum of the physical areas of the two nozzles must equal the 
nozzle area set by the engine control. 

The standard isentropic compressible flow equations were assumed to describe the 
flows through the two exhaust nozzles. Because of the assumption of an instantaneous 
pressure balance at the entrance to the augmentor (eq. (40)), it was possible to compute 
the low-pressure-turbine flow (eqs. (33) and (34)) without having to depend on volume 
dynamics downstream of the low-pressure turbine or resort to iterative procedures. 

The core augmentor equations (eqs. (41) and (44)) define the conditions at the entrance 
to the core exhaust nozzle and allow a direct calculation of the core nozzle area. 

In the nonaugmented mode, the physical exhaust nozzle area is constant. In order 
to provide a degree of freedom in matching digital nonaugmented data, variable flow co- 
efficients C d g and C d x g were assumed for the core and duct nozzles, respectively. 
Each flow coefficient was calculated from a linear function of its corresponding nozzle 
pressure ratio. The slope and intercept of each flow coefficient were adjusted to best 
match the digital data. 

During the augmented mode, the physical nozzle area was controlled so as to main- 
tain engine variables near their military power values. Therefore, the separate nozzle 
pressure ratios (and hence, flow coefficients) were nearly constant over the range of 
augmentation. In order to provide a degree of freedom in matching digital data during 
augmentation, a total nozzle flow coefficient C d N was calculated from a linear func- 
tion of the physical nozzle area. The slope and intercept were adjusted to best match 
the digital data. The following equations describe the TF30-P-3 exhaust nozzle model: 

< P / P >CN = ~ (49) 

*7 
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(50) 


A 8 = 


K Cn( w 4.1 + w BLLT + w f7 )V^ 


P 7 FN7 C d; 8 


r.- 


FN7 = 


0. 2588 if (P/P) CN < 0. 53 


r , -10.7143 ./ r . n 0. 2857 

McN V 1 - [ P/P JCN 


if (P/P) CN > o. 53 


(51) 


C d 8 = 1. 3625 - 0. 7 158(P/P) 


where 0. 825 < C d g < 1. 0 


A l. 8 " C d, N A N " A 8 where C d, N A N “ A N, MIN 


C , „ = 1. 049 - 1. 622X10" 4 
a, N JN 
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_ C d, 1. 8 A 1. 8 P 1. 7 FN17 
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dnV t i. 
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0.2588 if (P/P) DN — 0. 53 
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r / -10.7143 f r ->0. 2857 
[ P/P ] D N V 1 - [ P/P ]dn 




> (57) 


if (P/P) DN > 0. 53 


C d, 1. 8 = L 575 ' L 00(P/P) DN where °- 825 - C d, 1. 8 - L 0 


(58) 


Engine Dynamics 

The equations previously discussed describe the steady-state operation of the vari- 
ous engine components. However, a transient engine simulation must also consider the 
effects of fluid compressibility, fluid momentum, energy storage, and rotor inertias. 
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In the TF30-P-3 hybrid simulation, inter component volumes were assumed at en- 
gine locations where either (1) gas dynamics were considered to be important or (2) gas 
dynamics were required to eliminate the need for iterative solutions (fig. 3). Modified 
forms of the continuity, energy, and state equations (ref. 11) were written for each vol- 
ume and solved for the stored mass, temperature, and pressure in the volume. When 
mixing of gases was not involved, a simple first-order lag form of the energy equation 
was used. The time constant for each lag was estimated to be equal to the ratio of the 
stored mass in the volume to the flow through the volume. 

Because of the real-time requirement and the computing equipment limitations, gas 
properties such as specific heats were assumed to be constant over the range of engine 
operating temperatures and fuel-air ratios. Gas properties were adjusted at the engine 
design point (sea-level, static, standard-day, military power setting) to achieve the de- 
sired steady-state cycle balance. The following equations describe the intercomponent 
volume dynamics in the TF30-P-3 hybrid simulation: 
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The effects of fluid momentum on the TF30-P-3 engine dynamics were also included 
in the hybrid simulation. While these effects are primarily high frequency in nature, 
their inclusion in the simulation also serves to eliminate algebraic or iterative loops. 
The momentum equation was solved for the flow through (1) the combined fan core and 
low-pressure compressor, (2) the high-pressure compressor, (3) the combustor, (4) the 
fan duct, and (5) the duct augmentor. The following equations describe the TF30-P-3 
flow dynamics: 


1 / FID+LC J ° 


L ( P 2. 2’ " P 2. 2^ dt + w 2. 1, i 


3’ ‘ 


2. 2, i 


16 



'» - PM / t(P 4'- P 4 )dt + w 3,i 
\ / B 


( 86 ) 


w 


/AgA /*t v 

1,3 = H / ( p l.6'- p l.6 )dt+w 1.3,l 

\ L /D U 


(87) 


/AgA /-t . 

w 1.6 = (— ) J (P 1.7’- P 1.7 )dt + W 1.6,i 

V " / DAB U 


( 88 ) 


Probably the most significant factors in governing the transient behavior of a turbo - 
fan engine are the rotor moments of inertia. The rotor speeds were computed by using 
the conservation of angular momentum. The following equations are used in the 
TF30-P-3 hybrid simulation to compute the rotor speeds: 
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HYBRID SIMULATION 


The equations describing the TF30-P-3 mathematical model were implemented by 
using the Lewis Research Center's Analog and Hybrid Computing Facility. The Lewis 
EAI Model 690 Hybrid Computing System consists of an EAI 640 Digital Computer, an 
EAI 693 Hybrid Interface Unit, and an EAI 680 Analog Computer. Two EAI 231-R Analog 
Computers were also used in this simulation. 

The bulk of the calculations were performed on an EAI 680 and two 231-R Analog 
Computers. The analog computers perform all the operations characteristic of analog 
machines (i. e. , summing, integration with respect to time, multiplication, attenuation, 
univariate function generation, etc. ). The use of peripheral equipment such as x-y 
plotters and strip-chart recorders allows continuous monitoring by the user of computed 
variables. All the required analog computers can be tied together by means of a Cen- 
tralized Trunk System to allow the transmittal of information between analog consoles. 

The 640 Digital Computer was used primarily to perform the bivariate function gen- 
eration associated with fan, compressor, and turbine performance. The digital com- 
puter also provided teletypewriter output listings of selected engine variables. In order 
to minimize the core storage requirements and digital cycle time, scaled-fraction vari- 
a bles and arithmetic routines were used in the digital program. Scaled-fraction vari- 
ables are limited to values between ±1. The use of scaled fractions necessitates the 
scaling of digital variables in the same manner as the analog variables are scaled. 

Figure 8 illustrates the split of the computational load between the various com- 
puters used in the TF30-P-3 real-time simulation. The split of the computational load 
was based on the computing equipment complement of each console and a desire to mini- 
mize trunking. The control inputs to the TF30-P-3 simulation (i. e., w p4 , w p7 , w pl 
XBL7, XBL12, and A^) were provided by an independent digital computer representa- 
tion of the TF30-P-3 hydromechanical control system. That control simulation was im- 
plemented on the Lewis SEL 810B computer and is described in a separate report 
(ref. 10). 


Digital Program 

MAIN digital program . - The MAIN digital program performs the following tasks in 
the TF30-P-3 hybrid simulation: (1) reads unsealed component performance data and 
map scale -factor data by means of paper tape, performs the required map scaling, and 
stores the scaled map data in the MAP2 COMMON block which is shared with the bivari- 
ate function interpolation routine; (2) computes fan inlet conditions P 2 and T 2 from 

equations (1) to (5) and sets potentiometers on the analog to values corresponding to P 0 , 
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T and Pq; (3) accepts map input data from the analog by means of the analog-to- 
digital converters (ADC's); (4) performs the required function generation by means of 
calls to the bivariate function interpolation routine MAP2 (discussed in a separate sec- 
tion); (5) outputs map values to the analog by means of the digital-to-analog converters 
(DAC's); (6) provides teletype listing (if desired) of selected engine variables in both 
engineering and SI units. Additional ADC's are used to input analog data when obtaining 
teletype listings. 

Appendix B contains a FORTRAN listing of the MAIN digital program and a detailed 
discussion of its operation. Appendix C contains a description of the required format for 
digital input data (component performance data). Samples of the digital data format are 
also provided. 

MAP2 function routine . - Because of the nature of the fan and compressor perform- 
ance maps (figs. 3, 4(a), and 5(a)), rectilinear interpolation cannot be used to generate 
these functions. Therefore, a radial -interpolation, bivariate -function -generation rou- 
tine, MAP2, was developed at the Lewis Research Center for generating this type of 
function. The MAP2 routine is based on the MAPFUN routine described in reference 12. 
The MAP2 routine uses scaled-fraction arithmetic routines to minimize core storage re- 
quirements and execution time. Because of its favorable qualities, MAP2 was also used 
to perform the bivariate function generation associated with the turbine performance, al- 
though rectilinear interpolation could have been used. 

Appendix D contains a FORTRAN listing of the MAP2-MAP2L function generation 
routine and a detailed discussion of its operation. Also included in appendix D is an as- 
sembly language listing of MAP2-MAP2L. The assembly language version was used in 
the real-time simulation of the TF30-P-3 engine to further reduce the execution time and 
core storage requirements. 

The digital program (with assembly language MAP2-MAP2L) required approximately 
10 000 words of core storage. The resulting digital sampling interval (with steady analog 
inputs) was 4.65 milliseconds. 


Analog Program 

The analog portion of the TF30-P-3 hybrid simulation performs all the required 
computations except the bivariate function generation previously discussed. The analog 
computation load is split between the 680 Analog Computer and two 231-R Analog Com- 
puters. The 680 Analog Computer also serves as the interface between the engine simu- 
lation and the SEL 810B digital control simulation. Approximately 185 amplifiers and 
85 multipliers were required to perform the analog calculations. 

Appendix E contains the analog patching diagrams for the TF30-P-3 hybrid simula- 
tion. For the most part, conventional analog techniques were used to implement the 


19 



engine model. Scale factors are not shown on the analog diagrams. Appendix E also 
contains a list of the scale factors and a list of the potentiometer settings at the design 
point. The following sections describe two features of the analog portion of the simula- 
tion. 

Digital-to-analog multiplier interface. - One of the factors in determining the size 
(number of consoles) of an analog simulation is the complement of multipliers on each 
machine. On the 680 Analog Computer, a total of 24 multiuse multipliers are available. 
These multipliers may be used to multiply, divide, square, or take the square root. As 
a supplement to the multipliers, the 680 also has six digital-to-analog multipliers 
(DAM’s) which allow a digitally supplied variable to be multiplied by an analog variable 
with the product appearing on the analog console. Each additional product is obtained at 
the cost of two DAC channels which may not be used as conventional DAC's. Appendix E 
contains a list of the products obtained with the DAM feature in the TF30-P-3 real-time 
simulation. 

Analog function generation . - In order to decrease the digital sampling interval, the 
univariate function generation associated with the low- and high-pressure-compressor 
bleeds (figs. 4(b) and 5(b)) was performed on the 680 Analog Computer. Two 10-segment, 
variable -diode function generators (VDFG's) were used. The map input variables 

^c^LC,M anc * (^ c )hc,m were computed on the analog from equations (16) and (22), 
respectively, and transmitted to the digital by ADC's. 

Two 10-segment, fixed-diode function generators (FDFG's) were used on the 231-R 
Analog Computers to calculate the core and duct nozzle compressible flow functions FN7 
and FN17. These functions (eqs. (51) and (57)) assume a constant specific-heat ratio y 
of 1. 4. 


RESULTS AND DISCUSSION 
Steady State 

The TF30-P-3 hybrid simulation with digital control was evaluated over a range of 
steady-state operating conditions. The purpose of the steady-state evaluation was to test 
the various curve fits and assumptions in the hybrid model. 

With the simulation operating at the sea-level, standard-day, static, military power 
point, excursions along the operating line were made by ramping the PLA signal to a 
desired power setting. After the transient settled out, a steady-state listing of selected 
variables was obtained at the teletype. Table I contains the teletype data for the sea- 
level, static operating line for power settings of 15° (idle). 70°, and 120° (full augmenta- 
tion). 
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Figure 9 contains plots of a number of the hybrid-computer -simulated engine vari- 
ables over the range of power settings. For comparison, the corresponding values ob- 
tained from the baseline digital program were also plotted. In order to better match the 
digital data at the very low power settings, the 12th-stage bleed opening (XBL12) was 
limited to 70 percent of its fully open value (table 1(a)). 

Figure 9 shows good agreement between hybrid and digital data at PLA's above 30 . 
Below 30°, discrepancies occurred in the sensitivities of rotor speeds and compressor 
discharge pressure to the power lever angle. For power lever angles below 30 , the 
hybrid-computer -simulated fuel flow exceeded the baseline digital value. This resulted 
in discrepancies in rotor speeds, temperatures, and so forth, for this range of power 
lever angles. However, because of the nature of the control design, errors in the engine 
variables can also cause errors in control outputs such as fuel flow. The observed dif- 
ferences in fuel flow, rotor speeds, and compressor discharge pressure were probably 
due to inaccuracies in the engine model at the low power settings. 

Since almost all the previously discussed assumptions were based on sea-level, 
static, operating-line data, it was decided to evaluate the steady-state performance of 
the hybrid simulation at a different flight condition. The steady- state data acquisition 
was repeated for a 6. 096-kilometer (20 000-ft), Mach 1. 2 flight condition. For this 
flight condition, the 12th-stage bleed was closed for the entire operating line, and its 
effect on the compressor performance need not be considered. The fan inlet conditions 
for this flight condition were computed by the digital program and were P 2 = H* 205 
N/cm 2 (16. 252 psia) and T 2 = 320. 25 K (576. 45° R). Table II and figure 10 contain the 
teletype listings and comparison plots for this flight condition. 

For this flight condition, excellent agreement between hybrid and digital data was 
observed at the low power settings. Some discrepancies were observed in the range of 
power lever angles between 45° and 69°. Some unexplained irregularities were observed 
in the digital values of P g 3 , T 4 , Wg, and N L in this range of power settings. In gen- 
eral, however, the same basic trends were observed in both the hybrid and digital data^ 

Other discrepancies were observed in the augmentation range of PLA between 110 
and 120°. In the digital model, the nozzle area reached its maximum value only at a 
PLA of 120°, while in the hybrid model the area was near its maximum at 110 . This 
contributes to the observed differences in augmentor fuel flow since these control outputs 
are correlated in the digital control. 

The general agreement of hybrid and digital data shown in figures 9 and 10 substan- 
tiates the assumptions of (1) constant combustor efficiency, (2) linear pressure ratio 
temperature ratio relations, and (3) pressure-sensitive flow coefficients for the core 
and duct nozzles. The agreement of hybrid and digital data also substantiates the fan and 
compressor models used in the real-time simulation. 

The observed discrepancies at low power settings can be compensated for by manip- 
ulation of the compressor bleeds (if applicable). Perhaps a more satisfactory solution to 
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the problem of matching steady-state data at low power settings is to implement variable 
specific heats Cp at the compressor discharge stations. Experience with similar simu- 
lations indicates that linear Cp-T relations, used only in the load torque calculations 
(eqs. (90) and (93)) have strong effects on the speed split at low power settings and can 
be "tuned" to match low-power-setting data. This problem is eliminated completely in 
the simulation described in reference 7, where all fluid properties are assumed to be 
variable. Unfortunately, the necessary computing equipment and computing time were 
not available to do this in this simulation. 

The use of digitally implemented bivariate functions in the hybrid simulation did re- 
sult in a better wide-range match of steady-state data than previously obtained by the 
authors in all-analog simulations. 


Transients 


The TF30-P 3 hybrid simulation with digital control was evaluated for power lever 
slams (rapid increases) and chops (rapid decreases) between idle and military power 
settings and also between military power and full-augmentation power settings. Trans- 
ients were run only for sea-level, standard-day, static conditions. Data obtained with 
the hybrid simulation were compared with the corresponding baseline digital data. The 
power lever angle signal was ramped from its initial value to its final value at a rate 
corresponding to 125 degrees per second. Selected analog variables were recorded by 
using strip-chart recorders. Details of the TF30-P-3 control simulation are contained 
in reference 13. 

The simulation was operated at the idle power setting (PLA = 15°). Table 1(a) lists 
steady-state values for selected engine variables at that operating point. The power 
lever angle was then ramped to a value (67°) slightly below the military power setting 
(69 ). The 2° margin in power lever angle was chosen so as not to inadvertently initiate 
augmentation. 

Figure 11 shows the response of selected engine variables to this throttle slam. 

Also shown are the corresponding baseline digital simulation variables. The power lever 
angle reached its final value in approximately 0. 4 second. Significant increases in the 
combustor fuel flow (commanded by the digital control) were not seen until about 0. 7 sec- 
ond into the transient. 

The digital simulation of the fuel control (ref. 13) includes a secondary fuel mani- 
fold fill-drain calculation. The filling of the manifold was initiated when the commanded 
fuel flow reached 0. 1260 kilogram per second, which is approximately equal to the sea- 
level, static, idle fuel flow. The manifold was filled when the commanded fuel flow 
reached 0. 1638 kilogram per second. When the manifold was filled, the commanded fuel 
flow reached the engine without the lag associated with the fill dynamics. Figure 11 in- 
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dicates that the manifold filling took 0. 7 second in the hybrid simulation and 1. 0 second 
in the digital simulation. This time difference was a result of the fuel flow differences 
caused by the discrepancies in the engine models at the low power settings. 

Comparison with digital data indicates that the hybrid simulation variables re- 
sponded more rapidly to the PLA slam. The differences in the T 4 , P g 3 , and N H re- 
sponses were caused by the differences in fuel flow supplied by the control. 

The 12th-stage bleed doors which are open at the idle condition are closed during 
the acceleration when the low-pressure-compressor discharge static pressure exceeds 
1 481 p + 6. 254 N/cm 2 . Figure 11 shows that the hybrid simulation bleed doors closed 
approximately 0. 6 second before the digital simulation bleed doors. This was due to the 

faster response of P g 2 2 ^^rid simulation. 

The nozzle area was at its maximum value at the static, idle condition. It was re- 
duced to its minimum value (fig. 11 (e)) by the digital control at a rate consistent with the 
baseline digital simulation. 

The hybrid simulation was operated at the PLA = 67° condition. A throttle slam 
was then run with the PLA signal ramped to 120° at the same 125-degree-per-second 
rate. Table 1(c) lists steady-state engine variables at the full-augmentation condition. 
Figure 12 shows the resulting response of selected engine variables to this throttle slam. 
Also shown are the corresponding baseline digital variables. The fuel flows to the five 
augmentor zones are metered as functions of PLA, Njj, Tg, and P g , 3 . The nozzle 
area is increased by the digital control so as to maintain a scheduled ratio of P g 3 /Pg 
during the augmentor operation. This maintains all engine variables upstream of the 
augmentor approximately constant at their military power values. 

Figure 12 indicates that the timings of the zone fuel flow initiations in the hybrid and 
digital simulations are in general agreement. Differences of 0. 1 to 0. 2 second were ob- 
served, except for zone 4. Zone 4 fuel flow came on about 0. 5 second early apparently 
because of an incorrect hold calculation in the simulated augmentor control. This hold 
was intended to represent the time required to fill the zone 4 manifold. 

Figure 12(h) shows agreement between the nozzle area responses. The same 0. 1- 
to 0. 2-second delay for the digital response was observed. The agreement in the area 
responses indicates that the augmentor inlet and high-pressure-compressor discharge 
pressure dynamics are approximately the same in the hybrid and digital simulations. 

Figures 13 and 14 contain the hybrid and baseline digital engine responses to throttle 
chops from PLA = 120° to PLA = 67° and from PLA = 67° to PLA = 15°, respec- 
tively. Figure 13 shows agreement between the hybrid and baseline digital results for 
full augmentation to military chop. The timings of zone fuel flow shutoffs and the mag- 
nitudes of the selected engine variables are in agreement. Figure 13(i) does show 
slightly less damping in the hybrid simulation’s response of Pg. 
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Figure 14 shows excellent agreement between the hybrid and baseline digital results 
for the military-to-idle chop. The response of fuel flow (and hence, P g 3 and T 4 ) 
agrees since the line fill problem that existed for the corresponding acceleration does 
not exist at the start of the deceleration. This results in excellent agreement in the 
opening of the 12th-stage bleed when P g 2 2 = 1 - 543 p 2 + 6 ' 495 N / cm2 - The response 
of the temperature T 4 to the throttle chops and slams was obtained by decreasing the 
gain of the (WT) 4 integrator (eq. (66)) by a factor of 25. This was necessary to match 
the baseline digital response and was an attempt to represent heat-transfer effects in the 
combustor. These effects were considered in the baseline digital program but were not 
included in the hybrid model. 


Frequency Responses 

A comparison of the hybrid and baseline digital engine dynamics (independent of the 
control) can be accomplished by opening the loops between the engine and the engine con- 
trol. This was done with the hybrid simulation at the sea-level, static, PLA = 50° 
operating point. A constant combustor fuel flow was supplied to the engine with the 

bleed closed and the nozzle area at its minimum value. A sinusoidal fuel flow 
signal was superimposed on the constant fuel flow signal, resulting in sinusoidal oscilla- 
tions in engine variables such as rotor speeds and pressures. The amplitude of the 
driving signal was ±10 percent of the constant input fuel flow. Both the driving signal 
and a selected output signal were input to an analog transfer function analyzer which 
calculated the amplitude ratio of and the phase angle between the two signals. 

Figures 15 to 17 contain the frequency responses of the high-pressure-rotor speed, 
the high-pressure-compressor discharge static pressure, and the core augmentor inlet 
pressure, respectively, to fuel flow perturbations. Also plotted are the results from a 
similar procedure using the baseline digital simulation and a digital transfer function 
analyzer. The digital transfer function analyzer could not perform adequately above 
1. 3 hertz because of low signal levels. The analog transfer function analyzer had pro- 
visions for increasing the signal levels at the higher frequencies, allowing analysis to 
10 hertz. 

Figure 15 shows general agreement between the hybrid and digital responses of N 
Since the low-frequency dynamics are caused primarily by rotor inertia, agreement in 
the dynamics indicates that the rotor dynamics have been adequately modeled in the hy- 
brid simulation. The observed difference in the steady-state gain (normalizing factor) 
would seem to indicate a discrepancy in the steady-state engine models. However, the 
steady-state digital data presented in figures 9(a) and (b) show a sensitivity 


24 



AN h AN h /APLA 
Aw f4 Aw f4 /APLA 

of 0. 992 rpm/(kg/hr), which is close to the value obtained from the hybrid simulation 
data. Therefore, errors in the digital analyzer results are suspected. 

Figure 16 shows the frequency response of the high-pressure-compressor discharge 
static pressure. Significant differences were observed in both the steady-state gain and 
the dynamics. The hybrid data show a break frequency around 2 hertz with higher order 
dynamics above 1 hertz. The higher frequency dynamics may be attributable to the in- 
clusion of combustor momentum effects in the combustor. The digital response is flat 
with increasing phase lag, characteristic of a dead time. A study of the transient re- 
sponse of P g 3 to the PLA chop from 67° to 15° (fig. 14) shows a lag rather than a 
dead time between P g 3 and w p4 . Again, errors in the digital analysis are suspected. 

Figure 17 compares the hybrid and digital responses of P g to the w p4 oscilla- 
tions. Agreement in the steady-state gain was observed. The dynamics were in general 
agreement in the low-frequency range. However, the hybrid data show a resonant con- 
dition at about 3 hertz. The digital data indicate a similar condition above 1 hertz, but 
the inability of the digital analyzer to function above 1. 3 hertz prevented a real compari- 
son in that frequency range. The observed resonant condition may be attributable to duct 
dynamics since the pressure P g is assumed to be proportional to the duct augmentor 
pressure P^ g . 


CONCLUDING REMARKS 

A real-time hybrid simulation of the TF30-P-3 turbofan engine has been developed. 
The simulation is primarily analog in nature but does use the digital portion of the hybrid 
computer to perform bivariate function generation associated with modeling the perform- 
ance of the engine's rotating components. The digital sampling interval is approximately 
4. 6 milliseconds and does allow real-time simulation without compensation for the phase 
lag caused by the digital sample -and-hold process. 

The hybrid simulation was based on the engine model used in a baseline digital sim- 
ulation provided by the engine manufacturer. Differences between the hybrid and digital 
models were dictated by the real-time requirement. The hybrid model assumes constant 
fluid properties and combustor efficiencies to minimize analog equipment requirements. 
The hybrid model also includes momentum effects in the compressors and combustors 
which allow the use of flow as an input to the compressor performance maps without the 
need for iteration. 

Results from the hybrid and digital simulations were compared. These comparisons 
indicated that the real-time hybrid simulation does adequately match the baseline digital 
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simulation for both steady-state and dynamic conditions. Discrepancies at low power 
settings did result in faster hybrid responses for accelerations from idle to military 
power. These discrepancies could be minimized by considering the effects of variable 
specific heats on the engine balance at low power settings. 

The real-time hybrid simulation provides a valuable analytical tool for research 
studies involving the TF30-P-3 turbofan engine. In particular, digital controls may be 
developed and tested on the simulation prior to actual engine testing. Transient data ob- 
tained with the hybrid simulation and a digital control are presented. These data indi- 
cate that the simulation does provide a realistic test vehicle for the control. In fact, the 
generation of the data presented did result in the discovery of a number of "bugs" in the 
control which might otherwise have gone undetected. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 8, 1974, 

501-24. 
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APPENDIX A 


SYMBOLS 
2 2 

A cross-sectional area, cm (in. ) 

C d nozzle flow coefficient 

Cp specific heat, j/kg-K (Btu/lbm-°R) 

dt time increment, sec 

FN7 core nozzle flow function 

FN17 duct nozzle flow function 

general function, i = 1, 10 

fp flow parameter, kg-K-cn-r/N-rpm-sec (lbm- R-in. /lbf-rpm-sec) 

p 

g c gravitational conversion factor, 100 cm-kg/N-sec (386. 3 lbm-in. /lbf- 

HVF heating value of fuel, J/kg (Btu/lbm) 

Ah enthalpy drop, J/kg (Btu/lbm) 

hp enthalpy drop parameter, J/kg-K^^-rpm (Btu/lbm-°R^^-rpm) 

2 2 

I polar moment of inertia, N- cm -sec (in. -lbf-sec ) 

J work conversion factor, 100 N-cm/J (9339. 6 in. -lbf/Btu) 

K constant (appropriate units) 

l length, cm (in. ) 

M Mach number 

m second function map index, m = 7 to 12 

N rotor speed, rpm 

n first function map index, m = 1 to 6 

p 

P total pressure, N/cm (psia) 

P/P pressure ratio 

PLA power lever angle, deg 

Q torque, N-cm (in. -lbf) 

R a gas constant of air, N-cm/kg-K (in. -lbf/lbm-°R) 

T total temperature, K (°R) 

T/T temperature ratio 


V 


3 •} 

volume, cm (in. J ) 


W stored mass, kg (lbm) 

w mass flow rate, kg/sec (lbm/sec) 

XBL7 fraction of 7th-stage bleed 
XBL12 fraction of 12th-stage bleed 
XI ADC channel I input, I = 0, 23 

XSC scale factor on x map input 

x x input to bivariate map (normalized by design-point value) 

1 u 

x. array of x inputs for j value of y input, j = 1, 8 

YI DAC channel I input, I = 0, 23 
YSC scale factor of y map input 

y y input to bivariate map (normalized by design-point value) 

ZSC(l) scale factor on first function of x and y 
ZSC(2) scale factor on second function of x and y 
z l first function of x and y (normalized by design-point value) 

z 2 second function of x and y (normalized by design-point value) 

z^ array of first function outputs for j th value of y, j = 1, 8 

z 2 j array of second function outputs for j th value of y, j = 1, 8 

y specific -heat ratio 

6 total pressure (normalized by sea-level value) 

Tj efficiency 

9 total temperature (normalized by standard-day value) 

r time constant, sec 

Subscripts: 

B combustor 

BLHT high-pressure-turbine bleed 
BLLT low-pressure-turbine bleed 
BL7 7th-stage bleed 

BL12 12th-stage bleed 

BWHT portion of high-pressure-turbine bleed doing work 
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BWLT 

CAB 

CN 

c 

D 

DN 

F 

FAN 

FID 

FID+LC 

FZ1 

FZ2 

FZ3 

FZ4 

FZ5 

FOD 

H 

HC 

HT 

I 

i 

L 

LC 

LT 

M 

MIN 

N 

PR6 

s 

0 


portion of low-pressure-turbine bleed doing work 
core augmentor 
core nozzle 

corrected to inlet conditions 
duct 

duct nozzle 
fuel 

total fan 

fan hub section 

fan hub and low-pressure compressor 

augmentor zone 1 

augmentor zone 2 

augmentor zone 3 

augmentor zone 4 

augmentor zone 5 

fan tip section 

high speed 

high-pressure compressor 
high-pressure turbine 
inlet 

initial condition 
low speed 

low-pressure compressor 
low-pressure turbine 
map 

minimum value 
total nozzle 

core-duct pressure ratio 

static 

ambient 
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1. 3 fan tip discharge 

1. 6 duct augmentor inlet 

1. 7 duct nozzle inlet 

1. 8 duct nozzle throat 

2 fan inlet 

2. 1 fan hub discharge 

2. 2 low-pressure-compressor discharge 

3 high-pressure -compressor discharge 

4 high-pressure-turbine inlet 

4. 1 low-pressure-turbine inlet 

5 low-pressure-turbine discharge 

6 core augmentor inlet 

7 core nozzle inlet 

8 core nozzle throat 
Superscript: 

( inlet to mixing station 

FORTRAN symbols - MAIN program: 

2 

AN total nozzle area, in. 

9 

ANSI total nozzle area, m 

ETAI inlet efficiency 

I integer 

IERR error flag for linkage routines 

LX array containing number of points per curve for each map pair 
J integer 

JX map scaling index 

JY array containing number of curves for each map pair 

JZ1 map scaling index 

JZ2 map scaling index 

KX array containing x out-of-range counts for each map pair 

KY array containing y out-of-range counts for each map pair 
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M 

MAP2 

MAP2L 

MM 

MO 

N 

NCV 

NTBL 

NX 

NY 

PADR 

PLA 

PS22 

PS22SI 

PS3 

PS3SI 

PVAL 

PO 

POSI 

P2 

P2Q0 

P2SI 

P3Q22 

QRBADS 

QSHYIN 

QSRUN 

QSSECF 

QSTDA 

QWBDAS 

QWFRL 


integer 

bivariate function (first function) 
bivariate function (second function) 
function relay address 
Mach number 
map index 

number of curves for map being read in 
number of data points for map being read in 
number of points per curve for map being read in 
number of curves for map being read in 
array of addresses of digitally set potentiometers 
power lever angle, deg 

low-pressure-compressor discharge static pressure, psia 

2 

low-pressure-compressor discharge static pressure, N/cm 
high-pressure-compressor discharge static pressure, psia 

2 

high-pressure-compressor discharge static pressure, N/cm 
array of settings for digitally set potentiometers 
ambient pressure, psia 

2 

ambient pressure, N/cm 
fan inlet pressure, psia 
inlet pressure ratio 

2 

fan inlet pressure, N/cm 
high-pressure-compressor pressure ratio 
linkage routine for reading ADC's 
linkage routine for console selection 
linkage routine for setting logic mode 
linkage routine for setting time scale 
linkage routine for transferring DAC contents 
linkage routine for loading DAC registers 
linkage routine for setting function relays 


QWPR 

linkage routine for setting potentiometers 


SENSW 

logical array for testing sense switches 


TO 

ambient temperature, °R 


TOSI 

ambient temperature, K 


T2 

fan inlet temperature, °R 


T2QO 

inlet temperature ratio 


T2SI 

fan inlet temperature, K 


T22 

low-pressure-compressor discharge temperature, 

°R 

T22SI 

low-pressure-compressor discharge temperature, 

K 

T3 

high-pressure-compressor discharge temperature, 

°R 

T3SI 

high-pressure-compressor discharge temperature, 

K 

T4 

high-pressure-turbine inlet temperature, °R 


T4SI 

high-pressure-turbine inlet temperature, K 


T6 

core augmentor inlet temperature, °R 


T6SI 

core augmentor inlet temperature, K 


VALS 

array containing unsealed, normalized map data 


VI 

potentiometer settings, I = 1, 3 


WA2 

fan airflow, lbm/sec 


WA2SI 

fan airflow, kg/sec 


WAR2 

fan corrected airflow, lbm/sec 


WAR22M 

low-pressure-compressor map corrected airflow, 

lbm/sec 

WF4 

combustor fuel flow, lbm/sec 


WF4SI 

combustor fuel flow, kg/sec 


WF7 

core augmentor fuel flow, lbm/sec 


WF7SI 

core augmentor fuel flow, kg/sec 


WF17 

duct augmentor fuel flow, lbm/ sec 


WF17SI 

duct augmentor fuel flow, kg/ sec 


WR22SI 

low-pressure-compressor corrected airflow, kg/sec 

X I 

ADC input variables, I = 0, 23 


XNH 

high-pressure-rotor speed, rpm 
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XNL low-pressure-rotor speed, rpm 

XNLR2 fan corrected speed, rpm 

XSC scale factor for x breakpoints on map being read in 

XVALS array containing scaled x breakpoints for each map pair 

YI array containing DAC initial conditions 

YI DAC output variables, I = 0, 23 

YSC scale factor for y breakpoints on map being read in 

YVALS array containing scaled y breakpoints for each map pair 

ZSC(l) scale factor for first function in map pair being read in 

ZSC(2) scale factor for second function in map pair being read in 

ZVALS array containing scaled map outputs for each map pair 

FORTRAN symbols - MAP 2 (MAP2L) function routine: 

I index for x breakpoints 

IX array of last search indices for x breakpoints 

J index for y curves 

JY array of last search indices for y curves 

KX array of x out-of-range counts 

KY array of y out-of-range counts 

MAP2 bivariate function (first function) 

MAP2L bivariate function (second function) 

N map index 

NX array of number of points per curve 

NXP number of points per curve for map being read 

NY array of number of curves per map 

NYC number of curves for map being read 

XFRAC fractional location of XIN between XLO and XHI 

XHI (I + l)*"* 1 breakpoint on inferred YIN curve 

XIN x input to map 

XLO I m breakpoint on inferred YIN curve 

XVALS array of scaled x breakpoints for each map pair 
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YIN 

YINCR 

YVALS 

ZL 

ZR 

ZVALS 


y input to map 

fractional location of YIN between and (J + l)^ 1 curves 

array of scaled y breakpoints for each map pair 

map output for XLO, YIN inputs 

map output for XHI, YIN inputs 

array of scaled map outputs for each map pair 
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APPENDIX B 


MAIN DIGITAL PROGRAM - DESCRIPTION AND FORTRAN LISTING 

Description 

The primary function of the MAIN digital program is to read the bivaiiate function 
input variables on the ADC channels, to call the function generation routine MAP2, and 
to transfer updated function output variables to the analog computer through the DAC 
channels. 

The MAIN digital program and the function generation routine MAP2 ai e capable of 
handling as many as six pairs of bivariate functions, with each function described by as 
many as eight curves and 10 points per curve. A COMMON block MAP 2 is shared by 
the MAIN digital program and MAP2 (see FORTRAN listings). This COMMON block 
contains arrays for function data XVALS, YVALS, and ZVALS; search indices K 
and JY; argument dimensions NX and NY; and out-of-range counters KX and KY. 

The MAIN digital program reads the unsealed, but normalized, component peifoim- 
ance data and associated scale factors and stores the data in the VALS array. Appen- 
dix C contains samples of the required format for reading the data. The data in the 
VALS array are divided by the appropriate scale factor; and the results are stored in 
the XVALS, YVALS, and ZVALS arrays. The search indices and out-of-range counters 
are initialized, and the argument dimensions are stored in the NX and NY arrays. The 
functions of the search indices and out-of-range counters are discussed in appendix D. 

After scaling and storing the component performance data, the MAIN digital pro- 
gram reads scaled-fraction DAC initial- condition data. These data are stored in the 
Yi array. By initializing the DAC's with specified values, the user can check the analog 
portion of the simulation prior to entering the digital computation loop. Appendix C 
contains a sample of the required format for the DAC initial condition data. Unused 
DAC's are initialized to zero. 

Table III lists the scaled variables and associated DAC channels which are trans- 
mitted to the EAI 680 Analog Computer by the MAIN digital program. Since the desired 
initial condition corresponds to the engine design point, the initial values of the scaled 
map outputs should be equal to the inverse of the map output scale factors. 

After initializing the DAC’s, the MAIN digital program reads, from the teletype, 
unsealed values for the variables P Q , T Q , and M Q . The program calculates analog 
potentiometer settings corresponding to the variables P 2 , T 2 , and P Q for the specified 
flight condition. The analog portion of the simulation is discussed in appendix E. The 
potentiometers to be set by the digital computer are specified as input data by means of 
paper tape. The format used for this data input is illustrated in appendix C. 
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The main computation loop in the MAIN digital program starts at statement 30. At 
this point, the required analog data are read from the appropriate ADC’s and stored in 
the locations specified by the arguments used in the call to the QRBADS linkage routine. 
These arguments indicate that a total of 10 ADC's are to be read, starting with ADCO. 
The contents of these ADC’s are to be stored sequentially in the locations reserved for 
the variables listed in the SCALED FRACTION declaration, starting with the vari- 
able X0. 

Table IV lists the ADC’s and the associated scaled input variables received from the 
EAI 680 Analog Computer. The scale factors associated with the map inputs are equal 
to the products of the design-point values and the XSC (or YSC) scale factors specified 
with the map data. For example, at the design point, ADC channel 5 receives a scaled 
signal equal to 1. 0/1. 9045, or 0. 52507. 

The MAIN digital program then calls the function generation routine MAP2 with the 
appropriate input variables. The calling statement is of the form: 

z x = MAP2(n, x, y) 

The scaled map output values returned by MAP 2 are transferred through the DAC chan- 
nels to the analog computer. For certain of the functions, another called statement 

z 2 = MAP2L(m) is used. The reason for using a different calling statement is discussed 
in appendix D. 

After the DAC variables in table III are transferred, a sense switch on the 640 con- 
trol console is tested to determine if additional data are to be input from the analog. If 
the sense switch is depressed (analog should be in "HOLD” mode), the ADC channels 
listed in table V are read. Before teletype listing, the variables to be listed are un- 
sealed to allow listing with both engineering and SI units. After teletype listing, the 
program returns to statement 30, where updated ADC's are read. 


FORTRAN Listing 


DIMENSION VALS (248) ,ZSC (2 ) , PADR (3) ,PVAL (3) 
C*****ADC VARIABLES 


SCALED FRACTION XI ,X1 ,X2,X3,X4,X5 ,X6,X7,X8,XP,X! I, XI 1 ,X!2 
1 X14,X15,X16,X17 f X18,Xl?,X2l,X21 ,X22,X23 
C*****DAC VARIABLES 


,X! 


3 , 


SCALED FRACTION Y4 ,Y5,Y6,Y7, YF.YP.Yl I.Y1 1 
C*****0THER VARIABLES 

SCALED FRACTION XVALS Cl 1,8,6) , YVALS (8,6) , ZVALS (10, 8,12), YI (24), 

1 MAP2,MAP2L,VI,V2,V3,SSQRT * 

COMMON/MAPS /XVALS, YVALS, ZVALS, IX(6) ,JY<6) ,NX(6) ,NY(6) ,KX(6) ,KYC6) 
LOGICAL SENSW 
REAL Ml 
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C*****INPUT REAL COMPONENT MAP DATA AND SCALE FACTORS 
TYPE I 

I FORMAT (/3X,54HPLACE DATA TAPE FOR MAPS NO. 1-8 IN HSPTR. THEN R-S- 
IR./-) 

PAUSE 

DO 1 I Nr 1 ,4 

READ(4,3) XSC , YSC ,ZSC ( I ) 

3 FORMAT ( 4F1 1 .4) 

IFCN-4) 4,5,5 

4 NCV=8 
GO TO 6 

5 NCVr 7 

6 NTBL- I 7*NC V 

READ (4, 7) (VALS(I) ,1=1 ,NCV) 

7 FORMAT (8F8.2) 

Jr NCV+1 

READ (4, 8) (VALS(I) , Ir J , NTBL ) 

8 FORMAT (8F8.4) 

C*****SCALE COMPONENT MAP DATA 
DO 9 J=1 , NC V 

DO 9 1=1 ,8 

YVALS CJ , N )r VALS (J)/YSC 

jx= ncv+ 2 *( cj-n*8+n-i 

XVALSCI.J, N)r VALS (JX> /XSC 

9 Z VALS (I,J,N)rVALS CJX+1 ) /ZSC ( 1 1 
lX(N)r 1 

JY(N)r 1 
NX C N ) r 8 
NY ( N ) r NC V 
KX(N)r g 
KY( N)r g 
Ig CONTINUE 
NCVr 8 

NTBL = NCV*31 
DO 13 Nr 5 , 6 

R EAD (4 ,3) XSC , YSC ,ZSC(1),ZSC(2) 

READ(4,ll)(VALS(I),Irl ,NCV) 

JrNCV+1 

READ (4,1! )(VALS(I),Ir J.NTBL) 

11 FORMAT (1 gF7 .4 ) 

DO 12 Jr 1 ,NCV 

YVALS (J, NIrVALS (J) /YSC 

DO 12 Ir 1 ,1 g 

JX= NC V+3 g* ( J- 1 ) + I 

XVALS ( I , J, N) r VALS (JX) /XSC 

JZ 1 = JX+ 1 • 

ZVALS(I,J,N)rVALS(JZl)/ZSC(l) 

JZ2r JX+2 g 

ZVALS ( I , J , N+2) = VALS(JZ2)/ZSC(2) 

12 CONTINUE 
IX ( N) r 1 
JY ( N)r 1 
NX ( N ) r 1 g 
NY ( N ) = NC V 

13 CONTINUE 
C*****INITIALIZE dacs 

TYPE 14 

14 FORMAT (/3X,64HPLACE DATA TAPE FOR DAC INITIAL CONDITIONS IN HSPTR. 
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I THEN R-S-R./) 

PAUSE 

READ(4,15)(YI(I),I=1 ,24) 

15 FORMAT (5S8) 

C***+*I NITIALIZE ANALOG CONSOLE 
GALL QSHYIN(IERR,680) 

CALL QSRUN(IERR) 

CALL QSSECFC I ERR ) 

DO 20 0 M=1 ,24 

MM=5*M- 1 

CALL QWFRL(MM, .FALSE. ,IERR) 

200 CONTINUE 

C***** TRANSFER INITIAL DAC VALUES 
CALL QWBDAS (YI , 0 ,24 , IERR) 

CALL QSTDA 

C*****SPECIFY OPERATING CONDITIONS FOR ENGINE 
TYPE 16 

16 FORMAT ( /3X,33HTYPE DESIRED VALUES FOR P0,T0,M0./> 

ACCEPT 1 7 ,P0 ,T0 ,M0 

17 F0RMAT(F7.3,Fg.3,F5.2) 

T2Q0=<1 .+ ,2*M 0*M 0 ) 

P2Q0 = T2Q0**3.5 
ETAI=1 .0 

I F(M0 ,LE . 1 . ) GO TO 18 
ETA 1= 1 .0- .0 75* (M 0-1 .0)**1 .35 

18 T2=T2Q0*T 0 

P2 = P2Q0*P0*ETAI 
P0SI=P0*. 68948 
T0SI = T0*. 55555 
P2SI=P2*. 68948 
T2SI = T2*. 55555 
TYPE 1 9 P2 T2 

19 FORMAT ( /3x)5HP2 = ,F7.3,7H, T 2 = ,F 8 . 3 /) 

Vl = P2/40. 

V2 = T2/1 000. 

V3=P0/40. 

****SET POTS WHICH ARE SENSITIVE TO OPERATING CONDITIONS 
TYPE 20 

20 FORMAT (/3X, 55 HPL ACE DATA TAPE FOR POT ADDRESSES IN HSPTR. THEN R-S 
1-R ./) 

PAUSE 

READ(4,22)(PADR(I),I=1 ,3) 

22 FORMAT (8A5) 

PVAL ( 1 )=V1 
PVAL(2)=V2 
PVAL (3)= V3 
DO 25 1=1,3 

DO 23 J= 1 ,3 

CALL QWPR (PADR(I),PVAL(I),IERR) 

IF (IERR ,EQ. 1 ) GO TO 25 

23 CONTINUE 

TYPE 24 , PADR ( I ) , IERR 

24 F0RMATC3X, A4 ,9H IERR = ,11/) 

PAUSE 

25 CONTINUE 
TYPE 26 

26 FOR MAT ( 3X ,55HSLA VE CONSOLES 4 AND 5 TO CONSOLE 1. MANUALLY GO TO I 
1C.) 


38 



TYPE 29 

29 FORMAT ( /3X , 44HPR OCEED TO DYNAMIC 
PAUSE 

C*****READ ADC VALUES AND GENERATE MAP 

30 CALL QRBADS(X« ,l,m,IERR) 

Y7rMAP2 (1 ,XI,X1 ) 

Y4 = MAP2 <2,X0,Xl) 

Y5=M AP2 (3,X2,X3) 

Y6=MAP2 (A ,X4,X5) 

32 Y8=MAP2 (5,X6,X7) 


PART OF PROGRAM BY R-S-R./) 
OUTPUTS 


Y9=MAP2L(7) 

34 Y10-MAP2(6,X8,X9) 

Y 1 1 =MAP2L (8) 

C*****TRANSFER UPDATED DAC VALUES TO ANALOG 


35 CALL QWBDAS (Y4 , 4 , 8 , IERR ) 

CALL QSTDA 

C»****0 UTPUT UNSCALED DATA AT TELETYPE IF DESIRED 
J F( . NOT ,SENSW( ! ) ) GO TO 31 
CALL QRBADS(Xl«,l I, 14, IERR) 

VF4 = X1 0 
WF4=WF4*3. 

VF4SI=WF4*. 45359 


VF7=X18 

WF7=WF7*6. 

VF7SI=WF7*. 45359 
WF1 7= Xl 5 
WF1 7= WF 1 7*8 . 

WF!7SI=VF1 7*. 45359 

AN=X13 

AN=AN*I 600. 

ANSI=AN*. 00064516 
PS22= XI 4 
PS22= PS22*1 51. 

PS22S 1= PS22* . 68948 
PS3= XI 1 
PS3=PS3*400. 

PS3SI=PS3*. 68948 
T22 = XI 6 
T22=T22*1500. 

T22SI=T22*. 55555 
T3=X1 7 
T3-T3*2 108. 

T3SI=T3*. 55555 

T4=X12 

T4=T4*4I««. 

T4SI=T4*. 55555 

T6=X19 

T6=T6*2000. 

T6SI-T6*.55555 

VAR2=Y7 

VAR2 = WAR2*25 0. 

WA2= (WAR2*P2*I .5497) /SORT (T2) 
VA2SI=WA2*. 45359 
XNLR2 = X1 

XNLR2=XNLR 2*20000. 
XNL=.04391*XNLR2*SQRT(T2) 
WAR22M=X4 
WAR22M=VAR22M*40. 



VR22SI=WAR22M*. 45359 

P3Q22=Y6 

P3Q22- P3Q22*5 . 

XNH=X2 I 

XNH=XNH*2f»f f . 

PLA=>23 

PLA=PLA*I5I. 

TYPE 59 

99 F0RMATU8X,36HTF3i-P3 SIMULATION STEADY-STATE DATA//) 

TYPE 1M,PDSI,PI 

I if ^ FORMAT (5X,9HPi = ,F7 .3 ,9X, 7HN/SQ CM,7X,2H( ,F7.3,9X,8HPSI A ) 
TYPE 1«I,TISI,T« 

1®1 FORMAT (5X,9HTI = ,F7.2,9X,7HK ,7X,2H( ,F7.2,9X,8HR ) 

TYPE l B2 # M I 

1*2 FORMAT (5X,9HMI r ,F7.4) 

TYPE l 03 ,P2S I , P2 

I 03^ FORMAT (5 X»9HP2 r , F7 .3 ,9X, 7HN/SQ CM,7X,2H( ,F7.3,9X,*HPSIA ) 

TYPE 104,T2SI,T2 

1 04^FORMAT(5X,9HT2 = ,F7.2,9X,7HK ,7X,2H< ,F7.2,9X,PHR ) 

TYPE ! 15 , PLA 

105 FORMAT (5X, 9HPLA = ,F7 .2 ,9X, 3HDEG) 

TYPE 106,WF4SI,VF4 

1 FORMAT (5X, 9HVF4 = ,F7.4,9X,7HKG/SEC ,7X,2H( ,F7.4,9X,*HLBM/SEC) 
TYPE 1 07,VF7SI,VF7 

I 07^FORMAT(5X,9HWF7 r ,F7.4,9X,7HKG/SEC ,7X,2H( ,D7. 4,9X,8HLBM/SEC) 
TYPE 1 08,WF1 7SI,WF17 

1 18^ FORMAT (5 X, 9HWFI 7 = ,F7.4,9X,7HKG/SEC ,7X,2H< ,F7.4,9X,8HLBM/SEC> 

TYPE 1 09,X2 1 

119 FORMAT (5X f 9HXBL7 = ,S7) 

TYPE 110, X22 

1 1 • FORMAT (5X»9HXBL12 = ,S7) 

TYPE III, ANSI, AN 

1 1 1 ^ F0RMAT(5X,9HAN : ,F7 .5 , 9X,7HSQ M ,7X,2H( ,F7. 1 ,9X ,8HSQ IN ) 

TYPE 1 12,XNH 

112 FORMAT <5X,9HXNH = ,F7 . 0 ,9X ,3RRPM ) 

TYPE 1 1 3 ,XNL 

113 FORMAT (5X,9HXNL = ,F7 . 0 ,9X,3HRPM) 

TYPE 1 14,PS22SI,PS22 

1H FORMAT (5X,9HPS22 r ,F7.3,9X,7HN/SQ CM,7X,2H( ,F7.3,9X,8HPSIA ) 

TYPE 1 15,PS3SI,PS3 

1^5 FORMAT (5X, 9HPS3 = , F7 .2 ,9X , 7HN/SQ CM,7X,2H( ,F7.2,9X,8HPSI A ) 

TYPE 1 1 6, WA2S I , WA2 

II® FORMAT(5X,9HWA2 = ,F7.2,9X,7HK G/SEC ,7X,2H< ,F7.2,9X,8HLBM/SEC) 
TYPE 117 T22SI T22 

117^F0RMAT(5X,9HT22 = ,F7.2,9X,7HX ,7X,2H( ,F7.2,9X,SHR ) 

TYPE I18,T3SI,T3 
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118 FORMAT (5X,9HT3 = ,F7.1,9X,7HK 

n 

TYPE 1 1 9 ,T4SI ,T4 

119 FORMAT (5X,9HT4 = ,F7.1,9X,7HK 

1) 

TYPE 12«,T6SI,T6 

121 FORMAT (5X»9HT6 = ,F7.l,9X,7HK 

1 ) 

TYPE 121 , WR22S I , WAR22M 

121 FORMAT (5X»9HWAR22M = ,F7 .3 ,?X, 7HK G/SEC ,7X,2H( ,F7.3,9X,8HLBM/SEC) 

1 ) 

TYPE 122 , P3Q22 

122 FORMAT (5X,9HP3Q22 = ,F7.4) 

GO TO 3« 

END 


,7X,2H( ,F7.1,9X,RHR ) 
,7X,2H( ,F7.1,9X,8HR ) 
,7X,2H( ,F7.1,9X,8HR ^ 
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APPENDIX C 


DIGITAL INPUT DATA 

In the MAIN digital program, each map or set of component performance data is 
stored with an identifying map number in the XVALS, YVALS, and ZVALS arrays. The 
map numbers are required when calling the MAP 2 routine to allow the location of the 
appropriate blocks of stored data. Table VI contains a list of the map numbers for the 
TF30-P-3 component performance data. 

A basic requirement for the MAP 2 data is that each curve of a given function must 
have the same number of breakpoints. However, different maps may have different 
numbers of curves or breakpoints per curve. In the TF30-P-3 simulation, the fan and 
compressor maps are represented by eight points per curve, while the turbine maps re- 
quire 10 points per curve. All maps have eight curves except the high-pressure turbine, 
which has seven curves. 

Tables VII to XII illustrate the format used to input the component performance 
data. Tables VII to X contain the data for maps 1 to 4, respectively. In the TF30-P-3 
simulation, all component performance data are normalized by the sea-level, static, 

PLA = 69° value of the corresponding variable. Although normalized, the input data 
may exceed unity. In order to use the scaled-fraction function routine MAP2, these data 
must be scaled in the MAIN digital program. Therefore, the first line of data for each 
map contains scale factors XSC, YSC, and ZSC(l) which are used to scale the already 
normalized data. The second line of data contains the normalized y values. The re- 
maining lines of data contain the x and pairs for each of the ascending values of y. 
For example, the third line of data contains the normalized values of the first four x 
and Zj pairs for the lowest value of y. The fourth line of data contains the normalized 
map outputs corresponding to the last four x and z^ for the lowest value of y. The 
x breakpoints can be different for each of the curves. 

Tables XI and XII contain the data for map pairs (5, 7) and (6, 8), respectively. For 
convenience, data for pairs of maps describing the performance of a particular compo- 
nent are combined. In order to use the MAP2L feature described in appendix D, break- 
points for both maps in the pair must be common. For each map pair, the first line of 
data contains scale factors XSC, YSC, ZSC(l), and ZSC(2). The last scale factor is 
used to scale the second map output z 2 _ The second line of data contains the normalized 
y values. The remaining lines of data contain the x, z^ and z 2 values for each of the 
ascending values of y. For the case of map pairs, the x, Zj, and z 2 data are sepa- 
rated with a whole line used for each. 

The MAIN digital program accepts (by means of paper tape) the scaled-fraction 
initial conditions for all 24 (numbered 0 to 23) of the DAC's. The first five lines of data 
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shown in table XIII illustrate the format used to input these data. The last line of data 
in table XIII illustrates the format used to input the addresses of analog potentiometers 
to be set by the MAIN digital program for specified flight conditions. 
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APPENDIX D 


MAP 2 FUNCTION GENERATION ROUTINE - DESCRIPTION AND FORTRAN 
AND ASSEMBLY LANGUAGE LISTINGS 

Description 

MAP2 is a function generation routine developed for the EAI 640 Digital Computer 
for handling a special type of function of two variables. This special type of function 
contains curves (constant y values) which are not defined over the same range of the 
x variable. Examples of such functions are fan and compressor performance maps for 
turbofan engines. Normal rectilinear interpolation cannot be used for these functions so 
a special interpolation scheme is used in MAP2. 

MAP 2 can accommodate six pairs of independent input variables and twelve function 
output variables. That is, in addition to the six independent bivariate functions, six addi- 
tional functions can be generated. However, these additional functions each must have 
arguments and argument breakpoints common to one of the six independent functions. 

This feature is useful in modeling fans, compressors, and turbines whose performance 
is usually described by pairs of bivariate functions. 

In order to reduce storage requirements and execution time, MAP 2 was designed to 
use scaled-fraction arithmetic and scaled-fraction data. The tables of function data, 
map input and output variables, and the routine MAP2 itself must be declared as scaled 
fractions in any program or subroutine which loads the function tables or calls MAP2. 

The arrays IX and JY containing the search indices are initialized to 1 in the calling 
program. These search indices are updated during each execution of MAP 2 and are 
used to speed up the table search for subsequent calls to MAP2. Details of the table 
search and interpolation scheme are discussed in reference 12. 

The arrays KX and KY containing the out-of-range counters for each function are 
initialized to zero in the calling program. These counters keep track of the number of 
consecutive times an argument exceeds the range of the input data for its respective 
map. A limited number of out-of-range calls to MAP2 are allowed in order to permit 
short -duration increases in map inputs (caused by noise, etc. ) without problem interrup- 
tion. If the out-of-range count reaches 25 for either map input, the program is inter- 
rupted and subroutine MOOR (map out of range) is called. In the TF30-P-3 simulation, 
MOOR places the analog computer in "HOLD, " prints an error message at the teletype, 
and pauses to await the user's corrective action. The maximum allowable out-of-range 
count and the action taken by MOOR are arbitrary and may be modified for a particular 
simulation. 
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As previously discussed, six additional functions can be accommodated by MAP2 if 
they have arguments and argument breakpoints common to one ot the original six func- 
tions. This feature reduces the storage requirements for the input data arrays XVALS 
and YVALS. In addition, the execution time can be reduced by skipping the table search 
and going directly to the interpolation calculations. 

Both FORTRAN and assembly language versions of MAP2 have been developed. 
These may be used interchangeably. However, for the digital sampling interval lengths 
required for real-time engine simulation, the assembly language version must be used. 
Both versions have been written in such a way that overflows are not possible if the 
scaled-fraction data have been loaded correctly. 

In order to handle the additional functions for which the table search can be skipped, 
an additional entry point, MAP2L, was added to the assembly language version. The 
calling statement for these functions requires only the map number m. Since the FOR- 
TRAN compiler for the EAI 640 Digital Computer does not allow multiple-entry-point 
subprograms, a separate subprogram, MAP2L, was written. 


FORTRAN Listings 


C 


SCALED FRACTION FUNCTION M A P2 ( N , X I N , Y I N ) 

SCALED FRACTION X I N « Y I N , Yl NCR ,XH I ,XL0 , XFR AC , ZL , 7R , Y 1 « Y2 , 

1 MAP2.XVALS (1 Z,P,5) ,YVALS(P,S> ,ZVALS( 1 0,S, 12) 
C0MM0N/MAPS/XVALS,YVALS,ZVALS,lX(fi),JYC6) f NXC6) f NY(6),KXCS),KY(S) 

, .IF ENTRY IS FRON FAP2L LOCK UP OUTPUT ONLY. 

IFCN.GT.6) GO TO 43? 

NYC= NY ( N) 


NXP= Nxm 
I- IX c N ) 

J= JY ( N ) 

100 Yl = YI N-YVALS (J ,N) 

IF(Y1 . GT. . 00C00S) GO TO 110 
IF ( Y 1 . EQ . . 0 0 (1 3 0S ) C-O TO 120 
IFCJ.LE. 1 ) GO TO 143 


Jr J- 1 
GO TO 100 

1 1 0 Y2 = YI N-YVALS (J+l , N ) 

IF ( Y2. LT . . 0 2 3 3 ?S ) GO TO 1 ? 0 
IF(Y2.Ea.. 0C03C3) GO TO 130 


J= J+l 

IFCJ.GE.NYC) GO TO 150 


Yl- Y2 
GO TO 110 
120 YINCR= .P3000S 
GO TO 1 ? 3 
130 YI NCR = . S9999S 
GO TO I S' 3 
140 YI NCR: .00020S 
GO TO 160 
150 YI NCR: .??9?9S 
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J= J- 1 

160 IFCKYCN) .GE.25) GO TO 500 
KY (M)-KY (M)+l 
GO TO 200 

1S0 YI NCR: Y1 / (Y1-Y2) 

190 KY (M) = E 

20 0 XLO=XVALS( I, J, m+YlNCR*(XVALS(I , J+l ,N)-XVALS(I ,J,N>) 
IF(XI.'J.GT.XLO) GO TO 21 ? 

IF(XIN.EO.XLC) GO TO 220 
IF ( I .LE. 1 ) GO TO 24? 

1= 1-1 
GO TO 200 

21 0 XH I - XVALS ( 1 + 1 , J , N)+YI NCR* (XVALS ( 1+ 1 , J+ 1 , N) -XVALS ( I + 1 , J,N)) 
IF(XIN.LT.XHI) GO TO 2^0 
IFCXIN.EQ.XHI) GO TO 230 
I - 1 + 1 

IF(I.GE.NXP) GO TO 250 
XLO=XHI 
GO TO 21 0 
22 e XFRAC: .00000S 
GO TO 29? 

230 XFRACr . 99999S 
GO TO 29? 

240 XFRACn .30000S 
GO TO 260 
250- XFRAC: .99999S 
1 = 1-1 

260 Ir CKXCN). GE.25) GO TO 500 
KX(N)=KX( N)+l 
GO TO 30? 

2S0 XFRAC: CXI N-XLO)/(XHI-XLO; 

290 KX(N)=0 

3 0 0 I X C N ) = I 

JYC N) = J 

4 00 ZR = ZVALS ( I-rl , J,N)+YINCR*(ZVAL2(I+1 , J+ 1 , N ) - 7 . VALS C I + 1 , J,N) > 

ZL= ZVALS ( I , J , N)+YINCR*(ZVAL9( I , J+l , N ) -ZVALS ( I , J, N) ) 

M AP2= ZL+XFR AC* (ZR-ZL) 

RETURN 

500 CALL MOOR (N,xIN, YIN) 

RETURN 

END 


SCALED FRACTION FUNCTION MAP2LCN) 

SCALED FRACTION MAP2 ,MAP2L , XI N , YI N 

C XI N AND YIN ARE DUMMY ARGUMENTS AND DO NOT AFFECT RESULT. 

MAP2L=NAP2(N,XIN,YIN) 

RETURN 

END 
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SUBROUTINE MOOR (M,X!M,YIM) 

FCALEP FRACTION XTM,YIM 

C TEST FOR OPERATE MOPE 

CALL PRAMKILnC) 

■IFCILOC .FO. A) CALL OSH(IEPP) 

Type F XT M,vtm 

FOP NAT C /7HMAP no. ,13,20m IMPUTE OUT OF RAN0E/5HXIN - 
1 PH VIM = ,S7/) 

PAUSE 

IFCILOC .F0.4) CALL opnp ( TFRR? 

RETURN 

ENP 


Assembly Language Listing 


FCPP 1 : 
0 ^ P2 r 
00 P 03 : 
e<* 0 P 4 : 
P P P t*5 s 
AAPPS: 
P*»P P7s 


fl(* m<* : 





cox m on 

X A PS 


AMPS: 



0 0 7^0 

XVALS 

QCC 

4*(> 


Al (*: 



0 0 0^0 

YVALS 

qcc 

4« 


AMI 1 1 



PI 7M 

ZV*LS 

pep 

Of (» 


0 0 012 s 



0 0 0 0 S 

IX 

pec 

6 


00013s 



ppppf 

JY 

BSS 

6 


BM14: 



m(*$ 

MX 

BSS 

* 


0^015: 



pp p <*C 

MY 

BSS 

s 


0 ^ 0 1 ^ • 



mpps 

KX 

pec 

6 


00017s 




Kv 

BBS 

6 


0 0 0 1 s 

ppup a 


00000 


BEL 

P 


0001 ^: 





MAXE 

l^App f MAP 2 L 

0 0 02 0 s 





EXTERN 

MOOR 





* 

C OXPLET r 

FUNCTION 

GENERATION 

FMTPY POINT 


ppppp 

0 0 0 0 0 

0 0 0 0 0 

X AP2 

ADR 

<* 


0 0 023: 

F? AP! 

53777 

0 0 0 0 0 


LX 

*-l 

SET UP RETURN 

0 0 05>4 s 

{* (* n (*2 

POPPA 

0 0 0 0 ^ 


ICX 

A 


P M2>5: 

a p {* p 3 

51351 

?* (* 3 £ 4 


STX 

RETURN 



(* (* (* (*4 

141774 



LA 

*-4 

MAP MO. APR 

0fl?27: 

0 0 0 05 

2G44P 



SSM 


SET IMP. BIT 

0 0 02 «s 

{*(* r**S 

151334 

(* {*340 


STA 

TEX 1 



0 0 0 0 7 

57333 

P P?42 


LX f I 

TEX 1 


*003*: 

0 0 0 1 0 

51334 

{*<*344 


STX 

N 


00031? 

0 0 0 M 

71331 

0 0 342 


A OX 

TEX! 


PPF3P: 

p n f> jo 

1 45 33 0 

c* ?34£ 


LA, I 

TEX 1 


00033: 

0 0 0 1 3 

1 5133? 

00345 


STA 

XIX 


PPP5A ? 

{* O C» 1 A 

7 1 3?5 

n f*3A2 


A OX 

TEN ! 


00035: 


1 A53°5 

P f*34? 


LA, I 

TEN l 



{*?(* 1 c 

15133? 

p {*34f 


ST A 

YT M 


00037* 



* 

crj 

CONST A V TS 

A *’ R vvaLS 

APPRESS (X=N) 


00 03* 

00017 

1 4 0 v yy 

00747 


LA , x 

v X-l 


0003° 

0 ^ 00 ? 

161333 

0 ^353 


STA 

XC ■ 'T 


0 004 0 

00 no i 

1 4 0 y yy 

0*755 


LA , X 

V v_ J 


0 004 ! 

00 022 

I 6J33? 

0 0354 


STA 

YC NT 


0 0 04? 

n n no 3 

I 4°yyy 

0274 1 


LA f X 

my- ! 

CMRVFF/N A? 

00043 

0 0 no a 

161323 

0 0 3 4 7 


ST A 

mvc 

MyCr mv ( m ) 

0 p 044 

0 0 0?s 

| 42XXX 

0*733 


LA , X 

NY- I 

POIMTS/Ci'RVC 

PM^5 

PPP?6 

1 61 322 

0 0 35 0 


STA 

My p 

NXPrNX(M) 

0 0046 

0 0 027 

1 42YXX 

0*717 


LA,X 

IX- I 


^ pi 7 

0 0 0 3 0 

1*1321 

0 0351 


STO 

I 

IOLDr ty(mj 

0 004* 

0 0 03 1 

1 4?yyy 

02725 


L A ,X 

JY-J 

ACCUCr JCLP 

0 004° 

0 0 032 

22777 

0000! 


ncx 

1 

X= M- 1 

0 0 05 " 

0 0 033 

*65 0 0 



fx • 


A=M-1, X=JCLn 

0 0051 

P 0 P34 

31300 

0 0334 


M 

f* 1 A X *’ C V 


0 0 05* 

p P p 3 5 

26540 



FQ 



00053 

0 0 036 

1 6 1 3 0 3 

0034 1 


STA 

DEL 

:MAXNCV*(M- 1 ) 

00054 

0 0 03 7 

15 1273 

0 033? 


A 

YV 

00055 

0 0 040 

161302 

0034? 


STA 

TFPI 

r 

0 0 056 



* 

DETERYIMF 

Y IVTERVAL <X=Jr 

JCLD) 

00P57’ 

0 0 04 1 

141305 

00346 

LOOP1 

LA 

YIN 


0 005* 

0 0 04? 

177300 

0034 ? 


S,IX 

TEMI 

Y(J,N) 

0P05? 

0 n 043 

*7412 



snz 

« 

ACC=YI 

00060 

0 0 044 

41012 

0 0 056 


J 

L00P2 


0 0 06 1 < 

0 0 045 

*74 00 



sn 



0 0 0 $2 

00 046 

41 0^5 

001*3 


J 

STY I W C 

ACCr fl 

00063' 

0 0 047 

255 0 0 



FX 


ACC = J 

0 0064 

0 0 05 0 

121*67 

0 0 33 7 


C 

ONE 


00065' 

0 0 05 1 

*74 0* 



SG 



00066 

00050 

4 1 0*4 

00076 


J 

Y2SML 


00067! 

0 0 053 

*6S 0 0 



FX 


X- J 

00 063; 

00054 

2*777 

0 0 0 0 1 


ncx 

1 

Xr Jr J- 1 

00060: 

0 0 055 

4 1 764 

0 0 04 1 


j 

L00P1 

00070! 

0 0 05 6 

1 61265 

0 0343 

L00P2 

STA 

TFF2 

r Y1 

00 071 : 

0 0 057 

2*00! 

0000I 


ICX 

I 

X= J+l 

000^*! 

0 0 0 6 0 

141266 

0 0346 


LA 

YIN 

00073i 

0 0 0 6 | 

I 77261 

0 0 340 


S,!X 

TFf 1 1 

YCJ+I ,M) 

00074! 

0 0 06 ? 

27414 



SP7 


ACCrYP 

00075: 

0 0 063 

4 1 03 * 

00113 


J 

CALC I 

00076: 

00 064 

27404 



SP 


00077: 

0 0 0 65 

410 07 

0 0 074 


J 

YFPJ1 


0007*! 

00 066 

*65 0 0 



FX 


ACCr J+l 

0007 °! 

00067 

121260 

00347 


C 

NYC 

000*0: 

0 007 0 

0 7 4 0 4 



SL 


00 0*1 : 

0 0 071 

41010 

00101 


J 

Y2BIG 


0 0 0 ** • 

0007? 

2650 0 



EX 


Xr Jr J+l 

000*3; 

0 0 073 

41763 

00056 


J 

L00P2 

ACCr Y2 

000 * 4 ; 

00 074 

1 41 0 44 

0034 0 

YFQJ1 

LA 

NAXP’C 

r *77777 

p 0 p 35 « 

0 0 075 

41 0*5 

00 1*2 


J 

SETJ 


000*6; 

0007 c 

*6500 


Y2*YL 

FX 

X 

II 

000*7: 

00077 

*6740 



CLP 


0 0 0 ** • 

0 010 0 

41004 

00 104 


J 

TYC NT 


000*0; 

00101 

26500 


Y2BIG 

FX 

Xr J+l 

00 0O0; 

0 010? 

2*777 

0 0*0 01 


PCX 

1 

Xr J 

00001 ; 

00103 

1 4 I 035 

0 034 0 


LA 

t*AXI M C 

r f 77 77 7 

000O?: 

0 0104 

1 6!*51 

00355 

TYC MT 

STA 

VI NCR 


00003; 

00 105 

1 4 1 0 4 7 

* ^354 


LA 

YC NT 


00004 ; 

0 0106 

1*1*30 

0 0336 


C 

FAXCMT 


00005: 

. 00107 

074 04 



SL 
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p fl PPC 

^11? 

41154 

PP974 


J 

A^ORT 


** *Q7 

PP! i 1 

7 I ° 4 3 

P P 5 5 ^ 


AON 

VC NT 


ft (« (IQC7 

PPI 1° 

4 ] fl I 4 

p a 1 of: 


J 

FTORJ 


P (K POO 



* 

CALCULATE 

YIMCR 

(X=J+1, ACC 

- Y2) 

p P ! f t* 

pp ! 1 3 

9 » I fl J 


CALC 1 

TC A 


ACCr-Y2 

**1 ?! 

PP! !4 

1 5 1 °27 

PP3^3 


A 

TE1«2 

Yl 

*P1 *2 

PP! 15 

1-51225 

P P34R 


FT A 

TEP 1 1 

=Y!-Y2 

P * 1 p 3 

PPM £ 

2574 * 



CLR 



**1 *4 

P P 1 17 

9554 P 



ED 



**1 *5 

fn <* i 7 p 

1 4 1 223 

PP343 


LA 

TEr*2 


**| *5 

PP12 1 

I 1221 

P P3^? 


D 

TFMi 


**1 *7 

PPlf>? 

22777 

P P P P 1 

FETJ 

DCX 

1 

X=J 

**1 *P 

PP123 

1 <1232 

PP355 

STYIMC 

FT A 

YINCR 


PP! PQ 

PP!£M 

?574 * 



CLR 



PP1 1 P 

PP125 

1^1227 

P P3M 


FTA 

YCNT 

CLEAR COUNT 

PP 1 1 1 

^ ^ 1 ^ ^ 

51224 

PP35? 

FTOR J 

FTX 

J 


PP1 12 



* 

?ET UP XVALF AMD 

zvalf addreffef 

PPI 13 

P P 1 2 7 

22777 

PPPP1 


DCX 

1 

XtJ-1 

PP1 14 

pc* | J P 

?£5 p fl 



EX 



**115 

P*1 31 

15121* 

PP341 


A 

DEL 


0*115 

PP!3?. 

312*3 

PP335 



FAX'IPT 


PPI 17 

PP 133 

255 4 a 



ED 



PPI 1 « 

fl P 1 34 

1 512*5 

P P 342 


FTA 

TEM1 

: J,N 

PP! 1? 

PP135 

151174 

PP331 


A 

XV 


PP 12 P 

P F ! 3 6 

151222 

PP3^P 


FTA 

XI JN 

=X(0, J,N> 

PP121 

PP137 

151175 

PP335 


A 

f’AXMPT 


PP122 

p P 1 4 P 

151221 

p (* j $ i 


FTA 

XI JIM 

=XC*, J+l ,M) 

PP 123 

P P I 4 I 

1412*1 

P*34? 


LA 

TEM1 

- J, N 

0*124 

PPM? 

151 171 

PP333 


A 

TV 


0*125 

P P ! 43 

151217 

p P3K? 


FTA 

ZIJN 

-Z(*,J,N) 

0*125 

P P ! 4 4 

151171 

PP335 


A 

MAX H PT 


P P 1 2 7 

PPI 45 

I 51216 

PP363 


FTA 

ZIJ1N 

=Z(*,J+1 ,N) 

PP122 

PPI 4f 

532*3 

PP35! 


LX 

I 


p(*lpo 



* 

DFTFRU I ME 

X INTERVAL (X=I= 

I OLD) 

0P13P 

PPI 47 

147212 

PPJCl 

L00P3 

LA, IX 

XI JIM 

X(I, J+l , N ) 

PP131 

P P 1 5 P 

1 7721 * 

P P ? c p 


F , IX 

XIJN 

xn,j,N) 

PP132 

P P 1 5 1 

3 ! 2 M 

PP355 


M 

YIMCR 


P PI 33 

PP152 

1572*5 

PP3S? 


A, IX 

XI JN 

ACCrXLO 

PP 1 34 

P p 1 53 

121172 

PP345 


C 

XI N 


PPI 35 

? P 1 54 

27412 



SGE 



0*135 

P P 1 55 

41*12 

PPI £7 


J 

L00P4 


PP 1 37 

PP156 

274 *° 



FG 



P P 1 3 R 

PP157 

41*3! 

PF?I P 


J 

XLORH 

X=I 

PP 1 3° 

P P 1 P 

2.55* * 



EX 


ACCr I 

PP14P 

pmm 

121155 

P P337 


C 

ONE 


PP141 

F P 3 S? 

-774*2 



FG 



PP142 

PPI £3 

/i 1 *3* 

P P° 1 3 


J 

X2FPL 


PP143 

p P 1 £4 

255** 



EX 


X=I 

PP 1 44 

PPI 65 

22777 

PFPF! 


DCX 

1 

X=I=I-1 

PP 1 45 

p p I ££ 

41751 

F F 1 ^ 7 


J 

LOOP3 


PP14* 

P P 1 £ 7 

151157 

PP35^ 

L00P4 

FTA 

XLO 


PP147 

PPI 7P 

22**1 

PPPP1 


ICX 

1 

X=I+1 

PPI 4^ 

PPI 71 

1-4717* 

P P3K1 


LA, IX 

XI JIM 

XCI+1, J+l ,» 

(*<*14° 

PP 1 72 

1 771 55 

P P3^P 


F , IX 

XI JN 

XCI+1 ,J,M) 

PP15P 

PP173 

3 115? 

PP355 


r* 

YI NCR 


PPI51 

PPI 74 

157154 

PP3^P 


A, IX 

XI JN 

ACCrXHI 

PPI 52 

P P I 75 

i?i i5* 

P P345 


C 

XI M 


PP 1 53 

P PI 7^ 

2741 ^ 



fle 




PM54 

0*177 

M 051 

**73* 


J 

CALC 2 


00155 

(* Oio 3 p 

274 04 



OL 



00156 

ripo * i 

^1 007 

**21 * 


J 

XLORH 

X=I=I+1 

00157 

(1 1*0 $o 

°<5 0 0 



EX 


ACC- I-f -1 

00 1 5 <? 

* *9(*3 

121145 

**35* 


C 

MXP 


0 0 | So 

* *? *4 

274 04 



OL 



P0160 

* *o 05 

4101! 

* *21 £ 


J 

X2PI0 


P (* 1 ^1 


o cc 0 0 



FX 


X: 1= I+! 

PP162 

**2*7 

4 1760 

* * 1 £7 


J 

L00P4 

ACCrXHI 

001 *5 



3 fc 

XIN=XLO CX=n OR 

XIM=XHI 

(X=I=I+ 1 ) 

PP1 *A 

**?! * 

26740 


XIORH 

CLP 



PM *5 

**21 l 

161 146 

**357 


OTA 

XFRAC 


00J 6 C 

ffl * o 1 n 

41051 

* *243 


J 

RANOP 


PPJ67 

* *9 i 3 

PCS 0 0 


X20 M L 

FX 


X=I 

001 6 P 

**9 1 4 

2C740 



CLP 



001*0 

**215 

4 1 0 0/5 

** 22 ! 


J 

TXCMT 


00170 

* *9 1 £ 

2 C5 0 0 


X2B 10 

FX 


x= 1+1 

f «*1 71 

**7I 7 

22777 

****1 


OCX 

1 

X: I 

00 172 

**77* 

1 4 1 12 0 

**3^* 


LA 

MAXI NC 

= *77777 

PM 75 

* *99 ] 

161156 

**357 

TXCVT 

OTA 

XFRAC 


00174 


141151 

**353 


LA 

XCNT 


001 75 

**7£3 

121115 

**336 


C 

MAXCMT 


PP 1 76 

**224 

27404 



SL 



P PI 77 

**705 

41047 

**274 


J 

ABORT 


PP1 7P 

* *?o£ 

71 125 

**353 


AO!* 

XCNT 


PP1 7? 

**227 

41014 

**243 


J 

RAMDR 


PP1 «P 



* 

CALCULATE 

XFPAC 

(X=I+I , 

ACC=YHI) 

PM M 

**73 * 

171126 

* *356 

C*LC2 

<: 

XLO 


0 P 1 ^ 2 ! 

**231 

161 I 1 1 

* *3^2 


OTA 

TFM 1 

=XHI-XLO 

001*5; 

* *23? 

2*74 0 



CLP 



PP 1 *4 ; 

* *9 3 3 

2*54 0 



FQ 



PPI *5 ! 

**234 

141 I 1 I 

* *345 


LA 

XI M 


0 P 1 « 6 > 

* *2 35 

171121 

* * 35 6 


S 

XLO 


001*7: 

p*73£ 

1 1 1 04 

**342 


D 

TEM 1 


P01 **• 

* *23 7 

16! 120 

**357 


OTA 

XFRAC 


PP1 *0; 

* *24 * 

22777 

**** 1 


DCX 

1 

X=I 

PPIPP; 

* *74 I 

26740 



CLP 



00191 ; 

* * 2^2 

1 *1 11 1 

**353 


OTA 

XCNT 

CLEAR COUNT 

0P192: 



* 

OFT UP REENTRY CO^PITTOM*? 

(x=n 

001 95: 

* *24 3 

5110* 

**351 

RA NOR 

OTX 

I 


001°4: 

* *744 

2**00 



FX 


ACCHM = I 

P P 1 95 : 

**945 

55 077 

* 0344 


LX 

M 

X=MAP MO. 

00196: 

* *74 £ 

1 62XXX 

*971 7 


OTA,X 

IX- ! 

IX ( N ) - I 

«M Q 7: 

* *74 7 

1411 05 

**35? 


LA 

J 


PPIOC?. 

**75* 

1 62XXX 

*7725 


OTA ,X 

JY- 1 

JY ( N ) r J 

PP199: 

* *75 1 

14! 1 02 

* *353 


LA 

XC*JT 


0 02 0 P : 

* *757 

1 *2XXX 

*2747 


OTA ,X 

KX-! 

KX ( N ) r X C M T 

PP9M : 

**753 

141101 

**354 


LA 

YC mj 


P 02 02 : 

**754 

1 *2XXX 

*7755 


STA.X 

KY-1 

KY(M)=YCNT 

P 02 05 : 

**255 

55 074 

**351 


LX 

I 


PP2P4: 




CALCULATE 

FUNCTION OUTPUT 

(x- n 

PP2P5: 

**25 6 

1471 05 

**363 

INTRP 

LA, IX 

ZIJIM 

7 C I , J + 1 , m ) 

0 02 0 6 : 

* *25 7 

1 77! 05 

0 *3<T9 


o,ix 

ZI JM 

7(1, J,N) 

00207: 

** 2 <* 

51075 

**355 


M 

YI NCP 


0020 *: 

* *2 £ 1 

15710! 

**36? 


A, lx 

ZIJM 


0020 °: 

* *?££ 

1610C? 

* * 3 4 2 


OTA 

T£XI 

= 7L 

0021 0 : 

**?£3 

P 2 0 0 1 

* * * * 1 


ICX 

1 

Y- 1+ 1 

0 0 ° 1 1 : 

* *7£4 

1 47077 

0*3* 3 


LA, IX 

ZIJIM 

7(1+1 , J+l ,M) 
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**2 ]?: 

P09<5 

I 7707* 

0 *3*9 


S,IX 

ZIJN 

7(1+1 ,J,N) 

**213: 

0 09** 

31 0*7 

0*355 


p 

YI NCR 


**214: 

00957 

157073 

0*7*2 


A, IX 

ZIJN 

;7D 

**215: 

0 09 70 

1 71 0*2 

**342 


s 

TEP 1 

- ZP- 7 L 

**21 *: 

0 097 1 

31 *** 

**357 


P 

XFRAC 

‘ 

**217: 

0 02 72 

15 10*0 

0*34? 


A 

TEP I 

PAP OUTPUT 

**2! *: 

0 0273 

45 0 7 1 

0 03*4 


J.I 

RETURN 


**? 1 9 . 



* 

PAP OUT OF RANGE 



**22*: 

0 02 74 

*0XXX 


ABORT 

L 

POOR 


**22! : 

0 02 75 

yyyyy 

*0344 


APR 

N 


PP222 : 

0 0975 

xxxxx 

0 0345 


ADR 

XI M 


**2*3: 

0 02 77 

xxxxx 

00345 


ADR 

YI” 


**224: 

0 0 3 0 0 

0 0 0 0 0 

000*0 


ADR 

0 


**225: 

0 03 *1 

45 0*3 

0*3*4 


J,I 

RFTURN 

OUTPUT: ? 

**22*: 



* 

FUNCTION LOOKUP ONLY ENTRY 

POI NT 

**227: 



* 

FFT UP OUTPUT TABLE ADDRESSES 

**?2*: 

0 03 02’ 

0 0 0 00 

0 0 0 0 0 

PAP2L 

ADR 

0 


**22°: 

0 03 03 

53777 

00300 


LX 

*- ! 

SET UP RETURN 

**23*: 

0 0 3 *4 

22**1 

0 0 0 0 J 


icx 

I 


("0231 : 

0 03 *5 

51 057 

0*3*4 


STX 

RETURN 


0*232: 

0 03 0* 

141774 

0 03*2 


LA 

*— A 

PAP NO. APR 

00233: 

003 0 7 

2*440 



SSN 


SET INP. BIT 

00234: 

0 03 1 * 

1*1032 

0 0342 


STA 

TEP 1 


00235: 

0 031 1 

145031 

* *342 


LA, I 

TEP! 

ACCUPrM 

**23*: 

0 03 12 

1 71 0°5 

00337 


p 

ONE 

ACC: N- 1 

00237: 

0 03 1 3 

31021 

0 0334 


P 

PAXNCV 


0 023°: 

0 03 1 4 

9*540 



EP 



0023^: 

0 03 15 

1 *1 *°4 

0034 1 


STA 

DEL 


**24*; 

0 031 * 

141034 

0*35? 


LA 

J 


00241 : 

0 03 1 7 

1 7! 020 

**337 


s 

ONE 

ACC= J- 1 

00242 : 

0 0 32 0 

151021 

0*34 1 


A 

DEL 


002^3: 

0 0321 

31*14 

0*335 


M 

PAXNPT 


00244: 

0 0322 

2**4 0 



EQ 


ACC: J , N 

00245: 

0 0 323 

151010 

0 *333 


A 

ZV 


0024*: 

0 0324 

1*1 *3* 

0 03*2 


STA 

ZIJN 

:Z(0,J,M) 

00247: 

00325 

151*10 

00335 


A 

PAXNPT 


0P24«: 

0 032* 

1 *1 03* 

* 03*3 


STA 

ZIJ1N 

:7(0, J+l ,N) 

0 02 4° ! 

0 0327 

53 poo 

0*351 


LX 

I 


00250: 

0 033 * 

4172* 

0025* 


J 

IMTRP 


0 025! 



* 

COPPON REFERENCES 



R*«?52 

0*331 

XXXXX 

77777 

XV 

ADR 

XVALS-! 


00253: 

0 0332 

XX^XX 

0*737 

YV 

ADR 

YVALS-1 


00254 

0 0333 

XXXXX 

01017 

zv 

ADR 

ZVALS-1 


00255: 

0 0334 

0 0*10 


PAX MOV 

DEC 

a 


0025* 

0 0335 

0 0 0 1 2 


PAXNPT 

DFC 

10 


0*257 

0 033* 

0 0 03 1 


PAXCNT 

DEC 

25 


00252 

0 033 7 

0 0 0 0 1 


ONE 

OCT 

! 


0025° 

0 034 0 

11111 


PAX INC 

OCT 

77777 


0 02*0 

0 034 1 

C| 0 0 0 0 

0 0 0 0 0 

DEL 

P22 

I,* 


**2*1 

0 0342 

fl ffl fl 0 0 

0 0 0 0 0 

TE^l 

BBS 

1,0 


0 02 *2 

0 0 34 3 

0 0*00 

0 0 0 0 0 

TEP2 

BSS 

I.* 


**2*3 

0 0344 

0 0 0 0 * 

0 0 0 0 0 

N 

BSS 

1,V 


002*4 

0 0345 

0 0 0 0 * 

0 0 0 0 0 

XIN 

BBS 

1,0 


002*5 

0 034* 

0 0 ffl 0 0 

0 0 0 0 0 

YIN 

BBS 

l.p 


002** 

0034 7 

01 0 0 0 0 

0**00 

NYC 

BBS 

1,0 


002*7 

0 03*0 

0 0 0 0 * 

***** 

NXP 

BBS 

1 ,0 


002* Q 

0035 1 

0*0*0 

0 0 0 0 0 

I 

BBS 

1,0 


002*° 

00759 

* 0 0 * 0 

0 0 0 0 0 

J 

BBS 

1,0 
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(!P?7P: 

00353 

p 0 P p 1 

0 0 0 0 0 

XCNT 

BBS 

! # 0 

0f?71 ! 

ffl 

P 0 0 0 {i 

0 0 0 0 0 

YCM 

BBB 

1 , 0 

p <*? 72 : 

PP355 

0 0 0 *0 

0 0 0 0 0 

YIMCR 

BPS 

\ # 0 

PP273- 

P 035* 

0 0 0 p 0 

0 0 0 0 0 

XLf 

BBS 

1 # P 

PP2 7/ j 

*035 7 

p p a p p 

0 0 0 0 0 

XFRAC 

BSF 

! f 0 

PF?75: 

0 03** 

0 0 0 0 0 

0 0 0 0 0 

XIJM 

RSB 

I f P 

PP?7^: 

003*1 

0*000 

0 0 0 0 0 

XI J I N 

BFB 

1 , P 

PP277; 

PP3*? 

0 0 0 0 0 

0 0 0 0 0 

7IJ*> 

BBS 

I , P 

P P2 7^: 

0 03*3 

0 0 0 0 0 

0 0 0 0 0 

ZIJ1 M 

BBB 

1 t P 

(*f*?7R: 

003*4 

0 0 0 0 0 

00 0 0 0 

RETlfRM 

BPS 

1 # fl 

P PO C?(J , 



0 0 0 0 0 


END 

P 
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APPENDIX E 


ANALOG PROGRAM - PATCHING DIAGRAMS, POTENTIOMETER SETTINGS, 

AND SCALE FACTORS 

The bulk of the calculations in the real-time TF30-P-3 simulation are performed on 
one EAI 680 Analog Computer and two EAI 231-R Analog Computers. Figures 18 to 20 
contain the analog patching diagrams for those computers. 

Table XIV contains a list of potentiometer definitions and settings for the analog 
portion of the simulation. Table XV contains a list of scale factors for the analog vari- 
ables. 

The EAI 680 Analog Computer has a complement of 24 quarter -square multipliers. 
This complement can effectively be increased by six by using the DAM feature. Digi- 
tally generated variables may be multiplied, on the 680, by analog variables with the 
product available at the analog. Table XVI lists the DAM variables in the TF30-P-3 
hybrid simulation. Figure 18 illustrates the role of the DAM'S in the TF30-P-3 simu- 
lation. 
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TABLE I. - TELETYPE OUTPUT LISTINGS FOR SEA-LEVEL, 


STANDARD-DAY, STATIC OPERATING LINE - 
CLOSED- LOOP CONTROL 
(a) Power lever angle, 15° 


P 0 

= 

10.132 

N/SQ CM 

T 0 

= 

288.14 

K 

M 0 

= 

. 0 0 0 0 


P2 

= 

10.132 

N/SQ CM 

T2 

r 

288.1 4 

K 

PLA 

z 

15.06 

DEG 

WF4 

z 

.1260 

KG/SEC 

WF7 

z 

.0023 

KG/SEC 

WF1 7 

z 

.0 066 

KG/SEC 

XBL7 

z 

.00109 


XBL12 

r 

.68603 


AN 

r 

.65410 

SQ M 

XNH 

r 

9970. 

RPM 

XNL 

z 

3962. 

RPM 

PS22 

z 

12.915 

N/SQ CM 

PS3 

z 

32. 11 

N/SQ CM 

VA2 

z 

40.02 

KG/SEC 

T22 

z 

322.56 

K 

T3 

z 

425.2 

K 

T4 

z 

755.2 

K 

T 6 

z 

606.1 

K 

WAR22M 

z 

1 0,805 

KG /SEC 

P3Q22 

- 

2.3632 



( 

14.695 

PSIA ) 

( 

518.66 

R ) 

( 

14.695 

PSIA ) 

( 

518.66 

R ) 

( 

,2779 

LBM/SEC) 

( 

. 0051 

LBM/SEC) 

( 

.0146 

LBM/SEC) 


< 

1013.8 

SQ IN ) 

( 

18.731 

PSIA ) 

( 

46.58 

PSI A ) 

( 

88.24 

LBM/SEC) 

C 

580.62 

R ) 

C 

765.3 

R ) 

( 

1359.3 

R ) 

c 

1091.0 

R ) 

( 

23.823 

LBM/SEC) 


(b) Power lever angle, 70° 


P 0 

z 

1 0.132 

10 

z 

288.1 4 

M 0 

r 

. 0 00 0 

P2 

z 

10.132 

T2 

r 

288.1 4 

PLA 

r 

69.98 

WF4 

z 

.8421 

WF7 

z 

.5737 

VFI7 

r 

.0013 

XBL7 

r 

.00085 

XBL12 

r 

.00048 

AN 

r 

.38205 

XNH 

r 

14042. 

XNL 

r 

9336. 

PS22 

= 

53.124 

PS3 

r 

161 .39 

WA2 

r 

1 04.46 

T22 

r 

504.85 

13 

z 

695.3 

T 4 

z 

1276.0 

T6 

z 

807.2 

WAR22M 

z 

12.819 

P3Q22 

r 

2.8938 


N/SQ CM 
K 

N/SQ CM 
K 

DEG 

KG/SEC 

KG/SEC 

KG/SEC 


SQ M 

RPM 

RPM 

N/SQ CM 

N/SQ CM 

KG/SEC 

K 

K 

K 

K 

KG/SEC 


( 14.695 

( 51*. 66 

( 14.695 

C 518.66 

( 1.8566 

< 1,2648 

C .1029 


( 592.1 


( 77.051 

( 234.08 

( 230.30 

( 908.75 

( 1251.7 

( 2296.8 

( 1453.1 

( 28.261 


PSIA 7 
R ) 

PSIA ) 
R ) 

LBM/SEC) 

LBM/SEC) 

LBM/SEC) 


SQ IN ) 


PSIA ) 

PSIA ) 

LBM/SEC) 
R ) 

R ) 

R ) 

R ) 

LBM/SEC) 


TABLE I 


Concluded. TELETYPE OUTPUT LISTINGS FOR SEA-LEVEL. 


STANDARD-DAY. STATIC OPERATING LINE - 
CLOSED- LOOP CONTROL 


(<•) Power lever angle. 120° 


P0 

~ 

! fl.l 32 

N/SQ CM 

T0 

z 

28*. 1 4 

K 

M 0 

r 

.0000 


P2 

r 

10.132 

N/FQ CM 

T2 

- 

2*8 • 1 4 

K 

Pi. A 

Z 

1 19.95 

DEG 

VF4 

r 

.*403 

K G/FEC 

WF7 

r 

2.1331 

K G/SFC 

WF1 7 

r 

2.571* 

XG/SEC 

XRl. 7 

r 

. 00007 


XBLI2 

I 

. 00036 


AN 

T 

.65309 

FQ M 

XNH 

- 

14042. 

RPM 

XNL 

" 

033*. 

RPM 

PF22 

: 

53.124 

N/FQ CM 


r 

161.10 

N/FQ CM 

WAP 

- 

1 04.40 

KG/FFC 

T2P 

7 

*105*16 

K 

T X 

Z 

71 1.7 

K 

7 A 

- 

1 276. 0 

K 

Tf 

t 

*11,6 

K 

pp M 


1 2 . *23 

X G/FFC 

c I r. o p 


2 , **5 *> 



( 

1 4.605 

PFI A 

) 

( 

51*. 66 

R 

) 

( 

14.605 

PFIA 

> 

( 

51*. 66 

R 

) 

( 

1 .*526 

LBM/FFC) 

( 

4.746* 

LBM/FFC) 

C 

5 .6600 

LBM/FFC) 


C 

1012.3 

C Q IN 

) 

( 

77,050 

PFIA 

) 

c 

233.7* 

PFT A 

) 

( 

230.36 

LRM/FFC) 

( 

000 . 

P 

) 

( 

12 * 1.2 

P 

) 

c 

2206.* 

P 

y 

( 

1 460.0 

p 

) 

< 

2*. 271 

L q M/FFC) 


TABLE n. - TELETYPE OUTPUT LISTINGS FOR SIMULATED 
6. 096- KILOMETER (20 000-FT) ALTITUDE, MACH 1.2 
FLIGHT CONDITION - CLOSED- LOOP CONTROL 
(a) Power lever angle, 15° 


P0 

z 

4.660 

N/SQ CM 

Tl 

= 

24*. 64 

K 

Ml 

z 

1.1999 


P2 

z 

1 1.205 

N/SQ CM 

T2 

z 

320.25 

K 

PLA 

z 

15.01 

DEG 

WF4 

z 

.2470 

KG/SEC 

WF7 

z 

.0026 

KG/SEC 

WF17 

z 

.0013 

KG/SFC 

XBL7 

z 

-.00012 


XBL12 

r 

.00146 


AN 

z 

.3305! 

SQ M 

XNH 

z 

I ! 687 . 

RPM 

XNL 

r 

7029. 

RPM 

PS22 

r 

29.453 

N/SQ CM 

PS3 

r 

69.21 

N/SQ CM 

WA2 

z 

76.72 

KG/SEC 

T22 

T 

439.24 

K 

T 3 

z 

564.0 

K 

T 4 

z 

903.0 

K 

T 6 

z 

634.0 

K 

VAR22M 

z 

11.003 

KG/SEC 

P3Q22 

z 

2.2451 



( 

6.760 

PSIA ) 

( 

447.56 

R ) 

( 

16.252 

PSIA ) 

C 

576.45 

R ) 

( 

.5445 

LBM/SEC) 

( 

,005* 

LBM/SEC) 

( 

.0 029 

LBM/SEC) 


( 

512.3 

SQ IN 

) 

( 

42.71* 

PSIA 

) 

C 

1 00.39 

PSIA 

) 

( 

169.14 

LBM/SEC) 

( 

790.64 

R 

) 

( 

1015.3 

R 

) 

( 

1625.4 

R 

) 

( 

1141.3 

R 

) 

c 

24.257 

LBM/SEC) 



TABLE II. - Concluded. TELETYPE OUTPUT LISTINGS FOR SIMULATED 
6. 096 -KILOMETER (20 000-FT) ALTITUDE. MACH 1.2 


FLIGHT CONDITION - CLOSED- LOOP CONTROL 
(b) Power lever angle, 70° 


Pi 

- 

4*661 

N/SQ CM 

( 

6.761 

PSIA ) 

ti 

r 

248.64 

K 

( 

447.96 

R ) 

Mi 

z 

1 .1999 





P2 

z 

11.219 

N/SQ CM 

( 

16.292 

PSIA ) 

T2 

Z 

321.29 

K 

( 

976.49 

R 7 

PLA 

z 

69.98 

DEG 




VF4 

z' 

.7981 

KG/SEC 

( 

1.7996 

LBM/SEC) 

WF7 

z 

.991 1 

KG/SEC 

( 

1.2151 

LBM/SEC) 

WF 1 7 

Z 

-.1114 

KG /SEC 

( 

-.1109 

LBM/SEC) 

XBL7 

z 

.11197 





XBL12 

z 

.11173 





AN 

z 

.39264 

SQ M 

( 

608.5 

SQ IN ) 

XNH 

z 

14296. 

RPM 




XNL 

z 

9393. 

RPM 




PS22 

Z 

93.174 

N/SQ CM 

< 

76.977 

PSIA ) 

PS3 

z 

199.81 

N/SQ CM 

( 

225.97 

PSIA ) 

WA2 

Z 

119.21 

KG/SEC 

( 

231 .95 

LBM/SEC) 

T22 

Z 

939.34 

K 

( 

971.82 

R ) 

T3 

z 

763.8 

K 

( 

1375. f 

R ) 

T4 

Z 

1316.4 

K 

( 

2369.6 

R ) 

T6 

Z 

843.1 

K 

( 

1517.5 

R ) 

VAR22M 

Z 

12.991 

KG/SEC 

( 

27.758 

LBM/SEC) 

P3Q22 

z 

2.7816 








(c) Power lever angle 


120° 


PI 


4.661 

N/SQ CM 

( 

6.760 

PSIA ) 

TI 

r 

248.64 

K 

< 

447.56 

R ) 

Ml 


1.1999 





P2 

= 

11.219 

N/SQ CM 

( 

16.252 

PSIA ) 

T2 

• 

321.29 

K 

< 

576.45 

R ) 

PLA 

- 

119.99 

DEG 




VF4 

- 

.8111 

KG/SEC 

( 

1 .7662 

LBM/SEC) 

WF7 

- 

1 .8866 

KG/SEC 

( 

4.1594 

LBM/SEC) 

WF17 

• 

2.6294 

KG/SEC 

< 

5.7968 

LBM/SEC) 

XBL7 

- 

.11134 





XBL12 

z 

.11122 





AN 

z 

.69423 

SQ M 

( 

1114.0 

SQ IN ) 

XNH 

z 

14299. 

RPM 




XNL 

z 

9396. 

RPM 




PS22 

Z 

93.187 

N/SQ CM 

( 

76.995 

PSIA ) 

PS3 

Z 

199.73 

N/SQ CM 

( 

225.87 

PSIA ) 

WA2 

Z 

1 19.24 

KG/SEC 

( 

232.02 

LBM/SEC) 

T22 

z 

939.34 

K 

( 

970.82 

R ) 

T3 

Z 

763.8 

K 

C 

1375.0 

R ) 

T4 

z 

1319,6 

K 

< 

2368.1 

R ) 

T6 

Z 

843.1 

K 

c 

1517.5 

R ) 

WAR22M 

Z 

12.977 

KG/SEC 

{ 

27.729 

LBM/SEC) 

P3Q22 

z 

2.8193 







TABLE III. - DIGITAL-TO-ANALOG CHANNELS AND 


VARIABLES FOR TF30-P-3 


DAC 

channel 

Digital 

variable 

Engine variable 

4 

Y4 

Fan hub pressure ratio 

5 

Y5 

Low-pressure-compressor pressure ratio 

6 

Y6 

High-pressure-compressor pressure ratio 

7 

Y7 

Total fan airflow 

8 

Y8 

High-pressure-turbine flow parameter 

9 

Y9 

High-pressure-turbine enthalpy parameter 

10 

Y 10 

Low-pressure-turbine flow parameter 

11 

Yll 

Low-pressure-turbine enthalpy parameter 


TABLE IV. - ANALOG-TO-DIGITAL CHANNELS AND 
VARIABLES FOR TF30-P-3 


ADC 

channel 

Digital 

variable 

Engine variable 

0 

xo 

Fan tip pressure ratio 

1 

XI 

Corrected fan speed 

2 

X2 

Low-pressure-compressor corrected airflow 

3 

X3 

Bleed-shifted low-pressure-compressor corrected speed 

4 

X4 

Bleed-shifted high-pressure-compressor corrected speed 

5 

X5 

High-pressure-compressor corrected speed 

6 

X6 

High-pressure-turbine pressure ratio 

7 

X7 

High-pressure-turbine speed parameter 

8 

X8 

! Low-pressure -turbine pressure ratio 

9 

X9 

Low-pressure-turbine speed parameter 





TABLE V. - SUPPLEMENTAL ANALOG-TO-DIGITAL CHANNELS 


AND VARIABLES FOR TF30-P-3 


ADC 

channel 

Digital 

variable 

Engine variable 

10 

X10 

Combustor fuel flow 

11 

Xll 

High-pressure-compressor discharge static pressure 

12 

X 12 

High-pressure-turbine inlet temperature 

13 

X 13 

Nozzle area 

14 

X14 

Low-pressure-compressor discharge static pressure 

15 

X 15 

Duct augmentor fuel flow 

16 

X 16 

Low-pressure-compressor discharge temperature 

17 

X 17 

High-pressure-compressor discharge temperature 

18 

X 18 

Core augmentor fuel flow 

19 

X 19 

Core augmentor inlet temperature 

20 

X20 

High-pressure-compressor rotor speed 

21 

X21 

Fraction of 7th- stage bleed 

22 

X22 

Fraction of 12th- stage bleed 

23 

X23 

Power lever angle 


TABLE VI. - COMPONENT MAP DESIGNATIONS 
FOR TF30-P-3 


Map 

index, 

n 

Map output variable 

Equation 

1 

Fan corrected airflow 

6 

2 

Fan hub pressure ratio 

8 

3 

Low-pressure-compressor pressure ratio 

14 

4 

High-pressure-compressor pressure ratio 

20 

5 

High-pressure-turbine flow parameter 

28 

6 

Low -pressure-turbine flow parameter 

33 


Map 

Map output variable 

Equation 

pair 



index. 



m 



7 

High-pressure-turbine enthalpy parameter 

30 

8 

Low-pressure-turbine enthalpy parameter 

35 
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TABLE VII. - TF30-P-3 FAN CORRECTED AIRFLOW DATA 


1 .4934 

2 

.1342 

1 .0816 





. 32 #! 

.4268 

.5335 

.6402 

.7470 

.8537 

.9604 

1.0671 

.If #0 

.5840 

.3420 

.4134 

.4099 

.3770 

.4407 

.3587 

.4674 

.3339 

.4895 

.3171 

.5104 

.2964 

.5494 

.2379 

. 0900 

.6381 

.3992 

.4746 

.4368 

.4555 

.4714 

.4356 

.3024 

.4145 

.5309 

.3934 

.5724 

.3446 

.5994 

.2595 

.000 0 

.6597 

.4164 

.5589 

.5029 

.5221 

.5764 

.4810 

.6073 

.4590 

.6393 

.3966 

.6478 

.3642 

.6743 

.2920 

.0000 

.7354 

.4906 

.6437 

.5879 

.6078 

.6663 

.5619 

.7117 

.4983 

.7222 

.4611 

.7277 

.4237 

.7367 

.3158 

.0000 

.7917 

.5814 

.7471 

.6408 

.7328 

.7397 

.69 09 

.31 06 

.6325 

.8251 

.5905 

.8381 

.5061 

.8391 

.3634 

. 0000 

.8665 

.6858 

.8665 

.7592 

.8600 

.8236 

.8453 

.8791 

.8228 

.9630 

. 763 ! 

.9930 

.6718 

.9990 

.4326 

. 0000 

.9730 

.7792 

.9730 

.9485 

.9708 

1.0274 

.9643 

1 .0974 

.9505 

1.1503 

.9236 

I .2027 

.8354 

1.2187 

.9299 

.0000 1 

.0556 

.8516 

1.0556 

1.0314 

1.0543 

1 .2092 

1.0513 

1.2947 1 

.0465 

1 .3681 

1.0331 

1.4135 

.9993 

1.4310 

.6165 


TABLE VIII. - TF30-P-3 FAN HUB PRESSURE RATIO DATA 


1 .4984 

2 . 

1342 

1 .3599 

.3211 

.4268 

.5335 

.6402 

. 0 0 0 0 

.5135 

.3746 

.5135 

.4995 

.5135 

.5244 

.5131 

. 0 0 0 0 

.5444 

.3746 

.5444 

.4995 

.5444 

.5244 

.5435 

. 0 0 0 0 

.5893 

.4370 

.5938 

.5764 

.5933 

.6078 

.5929 

. 0 0 0 0 

.61 19 

.4495 

.6618 

.62 93 

.6 622 

.6663 

.6591 

. 0 0 0 0 

.621 0 

.5089 

.71 16 

.6928 

.7348 

.7807 

.7411 

. 0 0 0 0 

.6663 

.5269 

.7892 

.8236 

.83 72 

.8791 

.8417 

. 0 0 0 0 

.6799 

.6078 

.8585 

.9485 

.9551 

1 .0274 

.9718 

. 0 0 0 0 

.7252 

.6703 

.9410 

.1213 1 

.0806 

! .2192 

1 .1 074 


.7470 

.8537 

.9604 

1.0671 

.4370 

.5135 

.4745 

.5135 

.5494 

.5122 

.5744 

.5117 

.4370 

.5444 

.4745 

.5444 

.5494 

.5430 

.5994 

.5426 

. 5 ! 19 

.5961 

.5409 

.5938 

.6293 

.5911 

.6743 

.5870 

.5424 

.6668 

.5879 

.6654 

.6953 

.6563 

.7367 

.6504 

.5814 

.7234 

.6408 

.7302 

.8106 

.7425 

.8391 

.7434 

.6858 

.8204 

.7592 

.8304 

.9260 

.8454 

.9990 

.8522 

.7792 

.91 06 

.8646 

.9342 

. 0974 

.9832 

1.2187 

.9972 

.8516 

.9995 

1.0314 

1.0539 

.2947 

1.1323 

1.4310 

1.1650 


TABLE IX. - TF30-P-3 LOW -PRESSURE -COMPRESSOR 
PRESSURE RATIO DATA 


1.1915 

2 

.4246 

1 .8953 

.3637 

.6161 

.7274 

.8486 

.1747 

.3980 

.3136 

.3942 

. 40 ] ! 

.3824 

.4394 

.3592 

.2939 

.5496 

.4400 

.5344 

.5344 

.5 049 

.5784 

.4746 

.3971 

.6937 

.5608 

.6671 

.6603 

. 62 09 

.7066 

.5815 

.50 04 

.8263 

.6965 

.8165 

.8320 

.7789 

.8606 

.7088 

.6514 1 

.0235 

.8835 

1 .0223 

.9818 

.9192 

1 .0019 

.8244 

.7070 1 

.1068 

.9327 

1.0765 

1 .0082 

.9412 

1.0231 

.8442 

.7864 1 

.2130 

.9532 

1.0993 

1 .01 73 

.9495 

1.0314 

.8510 

.9262 1 

.3949 

1 .0096 

1.1614 

1 .0387 

.9692 

1 .0451 

.8620 


.9698 

. 994 ! 

1.0062 

1.0547 

.3336 

.3931 

.3786 

.3881 

.4559 

.3429 

.5640 

.1895 

.4648 

.5314 

.5115 

.5163 

.5956 

.4526 

.7 022 

.1895 

. 589 ! 

.6626 

.6378 

.6379 

.7254 

.5549 

.8579 

.1895 

.7333 

.8150 

.8069 

.8013 

.8698 

.6679 

.9214 

.1895 

.9197 

1.0113 

.9683 

.9571 

1.0041 

.7733 

1.0168 

.1895 

.9591 

1.0534 

.9966 

.9848 

1.0244 

.7896 

1.0327 

.1895 

.9748 

1.0701 

1.0069 

.9946 

1.0319 

.7952 

1.0406 

.1895 

1.0201 

1.1178 

1 .0336 

1 . 02 04 

1 .1468 

.8070 

1.0517 

.1895 



TABLE X. - TF30-P-3 HIGH-PRESSURE-COMPRESSOR 


PRESSURE RATIO DATA 


I .4352 

1 

.9045 

1 .7463 





. 66(56 

.7618 

.8570 

.9522 

1.0000 

1.0475 

1.0951 


.5 023 

.6636 

.5637 

.5612 

.5759 

.5361 

.5884 

.5075 

. 602 ! 

.4757 

.6093 

.4582 

.6211 

.4170 

. 63*4 

, 34<>2 

.5413 

.8033 

.6376 

.7052 

.6566 

.6790 

.6746 

. 6437 

. 691 ! 

.60 00 

.6986 

.5745 

.7097 

.5158 

. 73 02 

.3492 

.5920 

.9867 

.7277 

.8812 

.7542 

.8519 

.7790 

.8096 

.8009 

.7537 

.81 02 

.7209 

.8235 

.6437 

.8486 

.3492 

.6710 1 

.2679 

.8443 

1.1072 

,8777 

1 .0677 

.9049 

1 .0069 

.9268 

.9280 

.9358 

.8826 

.9476 

.7806 

.9724 

.3492 

.6925 1 

.3534 

.9128 

1 .2235 

.9534 

1.1861 

.9789 

1.1096 

.9979 1 

.0143 

1 . 0 050 

.9608 

1.0140 

.8445 

1.0406 

. 3492 

.7176 1 

.4495 

.9713 

1 .3381 

1.0172 

I .2993 

1.0409 

1 .2078 

1 . 0567 1 

.0967 

1.0621 

1 .0356 

1 .0696 

.9070 

1.0980 

.3492 

.7463 1 

.5368 

1 .0323 

1.4575 

1 .0747 

1.4009 

1 .0965 

1 .2965 

1.1059 1 

.1659 

1.1087 

1 .0967 

1.1120 

.9542 

1.1195 

.3492 


TABLE XL - TF30-P-3 HIGH- PRESSURE-TURBINE FLOW AND 


ENTHALPY PARAMETER DATA 


2 . 

150 

1 .7036 

1 . 

2828 

2.3539 




.8347 

.8688 

.9028 

.9369 

. 97 ! 0 

1 . 005 ! 

! .0391 

1.0732 



.2115 

.4230 

.6345 

.8460 

1.0575 

1 .2690 

1 .4805 

1 .6920 

1 .9035 

2 . 

1 .2155 

1 .2155 

1.2155 

.2155 

1 .2071 

1.1545 

1.0519 

.8851 

.5772 


1 .7372 

1 .6737 

1 .5536 

.370 0 

1.1137 

.8545 

.6214 

.4096 

.1977 


.2115 

.4230 

.6345 

.8460 

! . 0575 

I .2690 

1.4805 

1.6920 

1 .9035 

2 1 

1.1648 

1.1648 

1.1648 

.1648 

I .1571 

1.1096 

1 . 0070 

.8466 

.5580 


1 .6772 

1.6101 

1 .5007 

.3241 

1.0720 

.8192 

.6002 

.3954 

.1906 


.2115 

.4230 

.6345 

.8460 

1.0575 

1 .2690 

1 .4805 

1 .6920 

1 .9035 

2 .* 

1.1186 

I .1186 

1.1186 

.1186 

1.1103 

1.0590 

.9621 

.8017 

.5195 


1 .6242 

1 .5536 

1 .4477 

.2782 

1.0331 

.7909 

.5790 

.3742 

.1871 


.21 15 

.4230 

.6345 

.8460 

1.0575 

1 .2690 

1.4805 

1 .6920 

1 .9035 

2! 

1,0750 

1 . 075 0 

1 . 075 0 

. 0750 

1 .0667 

1.0198 

.91 72 

.7568 

.4874 


1 .5607 

I .5 042 

1 .3983 1 

1 .2358 

.9957 

.7556 

.5508 

.3629 

.1765 


.2115 

.4230 

.6345 

.8460 

1.0575 

1 .2690 

1.4805 

1 .6920 

1.9035 

2! 

1.0352 

1 .0352 

1,0352 1 

! .0352 

1.0262 

.9749 

.8723 

.7119 

.4490 


1 .5183 

! .4653 

1.3594 1 

1,1935 

.9604 

.7274 

.5296 

.3460 

.1694 


.2115 

.4230 

.6345 

,8460 

1.0575 

1 .2690 

1.4805 

1.6920 

1 .9035 

2! 

.9974 

.9974 

.9974 

.9974 

.9884 

.9364 

.8338 

.6735 

.41 05 


1 .4477 

1 .4 053 

1.3064 ] 

1.1475 

.9258 

.7062 

.5120 

.3354 

.1624 


.2115 

.4230 

.6345 

.8460 

1.0575 

1 .2690 

1.4805 

1.6920 

1.9035 

2! 

.962 7 

.9627 

.9627 

.9627 

.9538 

.8992 

.7953 

.6350 

.3848 


1 .3983 

1 .3594 

1.2688 1 

1.1122 

.8919 

.6709 

.4943 

.3177 

.1553 


.2115 

.4230 

.6345 

.8460 

1.0575 

1 .2690 

1.4805 

1 .6920 

1.0055 

2! 

.9300 

.930 0 

.9300 

,9300 

,921 1 

.8723 

.7697 

.6144 

.3656 


1 .3594 

1 .3241 

1.2358 ] 

1 .0805 

.8594 

.6497 

. 473 ! 

.3036 

.1483 



1 148 
Mlf 
MM 
I 1 43 
MM 
MM 
1 1 43 
MM 
MM 
I 143 

MM 

MM 

I I 43 

MM 
MM 
I 143 
MII 
000 0 
I ! 43 
0 0 0 0 
0000 
I I 43 
0 0 0 0 
0 0 0 0 


TABLE XII. - TF30-P-3 LOW- PRESSURE-TURBINE FLOW AND 
ENTHALPY PARAMETER DATA 


3 ,< 

6257 

I .85 I 8 A 

.6554 

3.2331 




• 3713 

.4629 

.5555 

.6481 

.8333 

.9259 

1.0185 

1.1111 



.3625 

.7251 

1.0877 

1 .4503 

I .8129 

2.1754 

2.5380 

2.9006 

3.2632 

3.6254 

2.6885 

2.6885 

2.6769 

2.6606 

2.6419 

2.6070 

2.5651 

2.5023 

2.3277 

. 0000 

2.6737 

2.4733 

1 . 972 ! 

1.5195 

1 .2340 

.9699 

.7436 

.501 1 

.2521 

. 0 0 0 0 

.3625 

.7251 

1 .0877 

1 .4503 

I .8129 

2.1754 

2.5380 

2.9006 

3.2632 

3.6254 

2 .2462 

2.2462 

2.2346 

2.2229 

2.2001 

2.1764 

2.1345 

2.0786 

1 .8389 

. 0000 

2.3441 

2 . §95 § 

1 .6553 

1 .3288 

1 .0705 

.8503 

.6304 

.41 06 

.1972 

. 0000 

.3625 

.7251 

1 .0877 

1 .4503 

1 .8129 

2.1754 

2.5380 

2.9006 

3.2632 

3.6254 

I .0621 

1 . 862 ! 

1 .8505 

1 .8342 

1 .0156 

1.7853 

1 .7458 

1 .6876 

1.3966 

. 0000 

2. §627 

1 .82 99 

1 .4775 

1.1962 

.9705 

.7339 

.5302 

.3427 

.1616 

. 0 0 00 

.3625 

.7251 

1 .0877 

1 .4503 

1 .8129 

2.1754 

2.5380 

2.9006 

3.2632 

3.6254 

1 .6061 

1 .5945 

1 .5828 

1 .5642 

1 . 539 ! 

! .5083 

1 .4618 

1.3966 

1 .1 630 

. 0 0 00 

1 .0687 

1 .6327 

! .3352 

1.0863 

.8713 

.6724 

.4784 

.2974 

.1325 

. 0 0 00 

.3625 

.7251 

1 .0877 

1 .4503 

1 .8129 

2.1754 

2.5380 

2.9006 

3,2632 

3.6254 

I .222 § 

1.2174 

1 .21 04 

1.1918 

1.1652 

1.1405 

! .0940 

1,0242 

,8379 

. 0000 

1 .5842 

1 .3643 

1.1154 

.8858 

.6883 

.5108 

.3556 

.2198 

,1 002 

. 0 0 0 0 

.3625 

. 725 ! 

1 .0877 

1 .4503 

I .8129 

2.1754 

2.5380 

2.9006 

3.2632 

3.6254 

1 . § 024 

1 .1824 

1 ,0707 

1.0591 

1.0363 

1 .0079 

.9683 

.9078 

.7448 

. 000 0 

1 .471 § 

1 .26 §9 

1 .0184 

.8050 

.6136 

.4526 

.3168 

.1939 

.0872 

. 0 0 00 

.3625 

.7251 

1 .0877 

1 ,4503 

1 .8129 

2.1754 

2.5380 

2.9006 

3.2632 

3.6254 

.9776 

.9776 

.9753 

.9590 

.9348 

.9124 

.8729 

.8147 

.6634 

. 0 0 00 

1 .341 7 

1.1380 

.9214 

.7177 

.5515 

, 404 ! 

.2780 

.1713 

.0775 

. 00 0 0 

.3625 

. 725 ! 

1 .0877 

1.4503 

I .8129 

2.1754 

2.5380 

2.9006 

3.2632 

3.6254 

.8961 

.8961 

.0845 

.8752 

.0519 

.8263 

.7867 

.7262 

.5935 

. 0 0 00 

I .2835 

1 . 070 ! 

.8406 

.6466 

.4995 

.3605 

.2554 

.1519 

.0614 

. 0 0 00 


TABLE XIII. - TF30-P-3 DIGITAL-TO-ANALOG 
CONVERTER INITIAL CONDITION AND 
POTENTIOMETER ADDRESS DATA 


.00000 

,00000 

.00000 

.00000 

.73550 

.52752 

.57264 

.92456 

.77955 

.42480 

♦21479 

.30936 

. 00000 

.00000 

. 0 00 00 

.00000 

.00000 

. 00000 

.0 0000 

, 0 00 00 

.00000 

.00000 

. 0 0 0 00 

.00000 


P 0 0 0 P 060 

PI 00 
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TABLE XIV. - POTENTIOMETER SETTINGS 


(a) EAI 680 Digital Computer 


Potenti - 

Definition 

Setting 

Potenti- 

Definition 

Setting 

ometer 



ometer 



POO 

P 2 40 

0. 3674 

P 19 

120/150 

0.8000 

P01 

2/3 

.6667 

P50 

1/100 

. 0100 

P02 

0. 7210/2 (eq. (10)) 

. 3605 

P51 

P s, 2. 2 /P 2. 2 

.9100 

P03 

0. 2796 (3/2) (cq. (10)) 

. 4194 

P52 

^15/518. 69 14.696/3 

. 8328 

P05 

1/3 

. 3333 

P55 

P LA./ 150 

. 4610 

P06 

P s,3 /3P 3 

. 3143 

P56 

PL W 150 

. 1000 

P07 

P s,3,i/ 100 

. 5774 

P60 

T 2 /1000 

. 5187 

P09 

w 2. l,i /20 ° 

.6112 

P61 

^518. 69/2000 

. 5092 

P10 

3(Ag c /O FID+IC /40 

. 8685 

P63 

V5 18. 69/1000 

.7202 

Pll 

2/3 

.6667 

P64 

5(0. 1989)/ 2 (eq. (17)) 

.4973 

P 12 

2/3 

. 6667 

P65 

(0. 84 - 0. 5 15)/ 2 (eq. (17)) 

. 1625 

P 13 

5(0. 06)/ 3 (eq. (19)) 

. 1000 

P67 

5(0. 4689 - 0. 19 89)/ 2 (eq. (17)) 

. 6750 

P14 

1/25 

. 0400 

P69 

1. 16/5 (eq. (17)) 

. 2320 

P 15 

*2.2,i /200 

.6112 

P70 

V 10 

.6555 

P 16 

3/40 

. 0750 

P71 

V 10 

.6555 

P 17 

(Ag c /i) uc /50 

. 1653 

P72 

T 3 ./1500 

. 6088 

P 19 

W F4, MAX 7 3 

. 5975 

P73 

3 °R A /V 3 

. 1506 

P20 

0. 1852/2 (eq. (23)) 

. 0926 

P74 

3/8 

.3750 

P21 

1. 0322/2 (eq. (23)) 

. 5161 

P79 

XBL7 . 

0 

P22 

5(0. 08032/2 (eq. (23)) 

. 2008 

P81 

VS 18. 69/1500 

. 5880 

P23 

5(0. 125 88)/ 2 (eq. (23)) 

. 3147 

P85 

0. 7 

.7000 

P25 

w f4 , i/3 

. 5975 

P88 

4(0. 055) (eq. (25)) 

.2200 

P26 

W F4, min 7 3 

.0948 

P99 

12(14. 696) ^/5/518. 69/ 27 

.6397 

P27 

w F4 /30 

0 

P100 

P 0 / 4 ° 

.3674 

P31 

0.7588/ 2 (eq. (12)) 

. 3794 

P 105 

1. 049/10 (eq. (54)) 

. 1094 

P32 

100 V518. 69/1000/6 (14.696) 

. 8184 

P 106 

A.. . /1600 

N. min' 

.3223 

P33 

3(0. 24 12)/ 2 (eq. (12)) 

. 3618 

P 107 

160(1. 622X10 -4 ) (eq. (25)) 

. 2060 

P35 

W F7, i /6 

0 

P 109 

XBL12 i 

0 

P36 

5/6 

. 8333 

Q02 

[( P P )lc/( P/ P ) LC. M-l 10 

. 1007 

P37 

5/6 

. 8333 

Q04 

C(P/ P)fip / ( P ' P) FID. m! 10 

. 1002 

P39 

0. 84/2 (eq. (17)) 

. 4200 

Q09 

a n i/1600 

.3223 

P40 

W F1. 7, i' 8 

0 

Q12 

l 2 (Ag c / 1 ) fjd+LC 400 

.3479 

P41 

5 8 

.6250 

Q14 

W 2.2.i 

.6216 

P42 

5/8 

.6250 

Q19 

3(a k c /Ohc /80 

.3106 

P43 

5, 8 

. 6250 

Q22 

(P P) HC /(P ' Pi MC. M 10 

. 1003 

P47 

15 

. 2000 

Q24 

PI. A 150 

.8333 



TABLE XIV. - Continued. POTENTIOMETER SETTINGS 


(b) EAI 231-4 Analog Computer 


Potenti- 

ometer 

Definition 

Setting 

Potenti- 

ometer 

Definition 

Setting 

P09 

1/2 

0. 5000 

Q12 

W l. 3, 1/ 300 

0. 3630 

P10 

( A gc/0 D / 375 

. 7224 

Q15 

3c p,2. 2'1.3 /10c p, 1.3 

. 4206 

Pll 

3/8 

. 3750 

Q16 

c p, 1. 3 1 ' 1. 3 /2c p, 1. 3 

.6980 

P12 

(Ag c /*) D /375 

. 7224 

Q20 

2(Ak c / ' /) DAB /75 

. 4014 

P13 

3/8 

. 3750 

Q21 

3/8 

. 3750 

P 14 

' 1. 3^ 2 

. 6980 

Q23 

1/3 

.3333 

P15 

15 V 2V l. 3 

. 1112 

Q24 

3c p, 1.3' 

.7177 

P 16 

c p, 3' 1. 3 //60c p, 1. 3 

. 0235 

Q26 

T l. 6 / 100 

.3732 

P19 

3/4 

. 7500 

Q39 

1/2 

. 5000 

P20 

2 ( A g c /0 DAB /75 

. 4014 

Q40 

8/150 

. 0533 

P21 

3/8 

. 3750 

Q41 

8l 'DAB HVF ■ !. 7' / ' 10 \ 1. 7 

. 3815 

P22 

w 1.6 ( i/ 300 

. 3630 

Q42 

3c p, 1. 6' 1. 7 /10c p, 1. 7 

. 3843 

P23 

1/5 

. 2000 

Q43 

1/6 

. 1667 

P24 

3c p, 2 

. 7176 

Q47 

4/5 

. 8000 

P25 

7 1.6 /100 

. 3732 

Q48 

0. 825(4/5) (eq. (58)) 

.6600 

P26 

10 V V 1.3 

. 1482 

P50 

15K D V 1. 3< A V V 4R A 

. 1592 

P27 

T 1.6.i /1000 

. 6640 

P51 

W 1.3.i/ 8 

. 3746 

P40 

1/5 

. 2000 

P52 

(WT)j 3 ./6000 

. 3315 

P41 

6}^ rj/50 

. 1537 

P54 

15K DAB V 1. fi (A "c / ' ,, DAB /4R A 

.0150 

P42 

1/5 

. 2000 

P55 

W 1.6.i/ 8 

. 3665 

P44 

125 V V 1. 7 

. 1029 

P60 

0.9yf2 K cn 

.6522 

P46 

4(1.575)/'50 (eq. (58)) 

. 1260 

P61 

135 4000 K dn /10 4 

. 4325 

P48 

0. 8(1.0) (eq. (58)) 

. 8000 

P62 

W. „ ./T50 

. 3386 



TABLE XIV. - Concluded. POTENTIOMETER SETTINGS 


(c) EAI 231-5 Analog Computer 





TABLE XV. - SCALE FACTORS 





TABLE XVI. - DIGITAL-TO-ANALOG 
MULTIPLIER UTILIZATION 
FOR TF30-P-3 


DAC channel 

DAM chan- 
nel (analog 
output) 

DAM variable 

Digital 

input a 

Analog 

input 

7 

19 

7 

P 2 (w c>FAN 

8 

20 

8 

n h^ w p'ht 

9 ! 

21 

9 


10 

22 

10 

n l ( Vlt 

11 

23 

11 

N L^ hp ^LT 


a See table III for list of digital input variables. 
"See figure 18 for source of analog input vari- 


ables. 


r Fan inside 
i diameter 



Figure 1. - Schematic of TF30-P-3 engine. 






















EAI 231 -R analog 
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Duct nozzle 
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Figure 8. - Split of TF30-P-3 simulation computation load. 






a □ □ 


O Hybrid simulation 

□ Base'ine r!iu«td I simulation 


fdi High -pressure-turbine inlet temperature as function of power lever a n cj t e . 


I I L 1 L I I 

iei High-pressure-compressor discharge static pressure as function of power 
lever angle, 


20 ' 30 40 30 ftO H\ 

Power lever angle, PLA, 1c-' 

if> Total fan airflow as function of power lever angle. 

Figure 9. - Continued. 
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Total nozzle area, A N , nr Duet augmentor fuel flow, w pl j, kg/sec Core augmentor fuel flow, w p7 , kg/sec 



(g) Core augmentor fuel flow as function of power lever angle. 



Ih) Duct augmentor fuel flow as function of power lever angle. 
.7 


.6 



Power lever angle, PLA, deg 
(i) Total nozzle area as function of power lever angle. 
Figure 9. - Concluded. 
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(g) Core augmentor fuel flow as function of power lever angle. 
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(h) Duct augmentor fuel flow as function of power lever angle. 
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(i) Total nozzle area as function of power lever angle. 


Figure 10. * Concluded. 
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(a i Power lever angle. (Full scale. 150°. > .0 Power lever angle. 'Full scale. 130° 



lei Total nozzle area. (Full scale. 1.032 m?. ' <j. Low-pressure-rotor speed, 'lull scale. 20 000 rpn. . 

Figure 14. - Comparison of hybrid and baseline digital simulation responses to throttle chop from power lever angle of 67° to 13 J at 123 degrees per second. 



O Hybrid simulation (normalized by 
1. 116 rpm/kg-hr" 1 ) 

Baseline digital simulation data (normalized 

by 0.7775 rpm/kg-hr m 



Frequency, Hz 

Figure 15. - Comparison of hybrid and baseline digital frequency responses of high - 
pressure-rotor speed to fuel flow oscillations. Sea-level, standard-day, static con 
ditions; power lever angle, 50°. 

O Hybrid simulation (normalized by 
0.0458 N-cm‘^/kg-hr‘1 

Baseline digital simulation (normalized by 

0.0341 N-cm"^/kg-hr"l) 





Frequency, Hz 

Figure 16. Comparison of hybrid and baseline digital frequency responses of high- 
pressure-conpressor discharge static pressure to fuel flow oscillations. Sea-level 
standard-day, static conditions; power lever angle, 50°. 



O Hybrid simulation {normalized by 
0. 00305 N-cm’^kg-hr" 1 ] 

Baseline digital simulation (normalized 

by 0.00311 N-cm' 2 /kg-hr"M 



.1 1 
Frequency, Hz 


gure 17. - Comparison of hybrid and baseline digital frequency responses of core 
augmentor inlet pressure to fuel flow oscillations. Sea-level, standard-day, 
static conditions; power lever angle, 50°. 
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(b) Low pressure compressor. 




PB4 
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Figure 19. - Patching diagrams for EAI 231 4 Analog Computer. 
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Figure 19. - Continued. 
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Figure 20. - Patching diagrams for EAI 231-5 Analog Computer. 



(b) High-pressure turbine. 
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